Superconductivity close to the Mott state: 
From condensed-matter systems to 
superfluidity in optical lattices 



Karyn Le Hur^ and T. Maurice Rice^ 

^Department of Physics, Yale University, New Haven, CT 06520, USA 
^ Theoretical Physics, ETH Zurich, CH-8093 Zurich, Switzerland 



Abstract 

Since the discovery of high-temperature superconductivity in 1986 by Bednorz and 
Miiller, great efforts have been devoted to finding out how and why it works. From 
the d-wave symmetry of the order parameter, the importance of antiferromagnetic 
fluctuations, and the presence of a mysterious pseudogap phase close to the Mott 
state, one can conclude that high-Tc superconductors are clearly distinguishable 
from the well-understood BCS superconductors. The d-wave superconducting state 
can be understood through a Gutz wilier- type projected BCS wave-function. In this 
review article, we revisit the Hubbard model at half-filling and focus on the emer- 
gence of exotic superconductivity with d-wave symmetry in the vicinity of the Mott 
state, starting from ladder systems and then studying the dimensional crossovers to 
higher dimensions. This allows to confirm that short-range antiferromagnetic fluc- 
tuations can mediate superconductivity with d-wave symmetry. Ladders are also 
nice prototype systems allowing to demonstrate the truncation of the Fermi surface 
and the emergence of a Resonating Valence Bond (RVB) state with preformed pairs 
in the vicinity of the Mott state. In two dimensions, a similar scenario emerges 
from renormalization group arguments. We also discuss theoretical predictions for 
the d-wave superconducting phase as well as the pseudogap phase, and address 
the crossover to the overdoped regime. Finally, cold atomic systems with tunable 
parameters also provide a complementary insight into this outstanding problem. 
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1 Introduction 



Shortly after the discovery of high temperature superconductivity in cuprates 
[1], Anderson [2] suggested that the key to this perplexing phenomenon resides 
in the large repulsive (positive) interactions in the copper oxide planes. As a 
matter of fact, at zero hole doping, the single-band Hubbard model captures 
the insulating behavior of the parent compounds. The origin of this insulating 
behavior was described many years ago by Nevill Mott [3] as a correlation 
effect. In the ground state, each Cu site would be in a (f configuration open- 
ing a charge gap of the order of the Hubbard interaction U. Essentially, the 
electrons become localized as a result of the Hubbard interaction, and their 
spins form an antiferromagnetic arrangement (the Neel state) as a result of 
the virtual hopping of the antiparallel spins from one Cu ion to the next (the 
parallel configuration being disallowed by the Pauli exclusion principle) [4] . It 
is relevant to observe that the layer symmetry breaks the degeneracy of the 
d orbitals down to a single d^-^-y-i orbital, so that orbital issues which cause 
complications in most other transition metal oxides do not play a role here. 
The energy gain due to this ordering is known as the superexchange energy 
J ~ At^ /U, where t is the hopping energy between Cu ions. 

The Hamiltonian of the Hubbard model takes the explicit form 

H^-tYl {dldjs + h.c.)+Uj24^di^4idii. (1) 

{i;3),s i 

Here, i and j run over lattice sites, (...) denotes nearest neighbors, and s the 
spin of an electron. The first part describes the hopping of electrons between 
adjacent sites on the lattice (formed by the Cu ions) and the second part the 
on-site repulsion U > Q oi electrons at the same site. One electron per lattice 
site corresponds to half-fiUing (undoped case) . The strong-coupling limit of the 
Hubbard model results in the t — J model. The two-dimensional (2D) doped 
Hubbard model has so far resisted a definitive solution, primarily because its 
spins and holes are highly entangled with no obvious small parameter to sep- 
arate them. Whether the 2D Hubbard model describes the physics of high-Tj. 
cuprates remains a great subject of interest. In the cuprates, upon doping 
with holes the antiferromagnetism is rapidly destroyed and above a certain 
level superconductivity irrefutably occurs with dx2_y2 pairing symmetry, as 
conclusively demonstrated by phase sensitive experiments [5,6]. In fact, the 
low-energy density of states can be directly accessed via scanning tunnel spec- 
troscopy experiments [7] . Such higher angular momentum pairs have led to the 
speculation that the pairing in the cuprates have a different origin from that of 
the well-understood Bardeen-Cooper-Schrieffer (BCS) superconductors [8]. To 
answer this question, one certainly needs to understand the normal phase at 
intermediate doping levels between the magnetic and d-wave superconducting 
phases (consult Fig. 1). 
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Fig. 1. Schematic phase diagram of high- Tc superconductors. 

There are experimental signs of a spin gap opening below a crossover tempera- 
ture T* [9,10]. The ultimate nature of the underlying quantum ground state in 
this portion of the phase diagram, commonly called the 'pseudo-gap' regime, is 
an intriguing puzzle. In 1987, Anderson [2] postulated a description based on 
a spin liquid of spin singlets which is lightly hole doped. Rather than forming 
a fixed array of dimers, strong quantum fluctuations lead to a superposition 
of singlet configurations, i.e., the valence bond singlets resonate between dif- 
ferent configurations. This is the Resonating Valence Bond (RVB) theory that 
predicts the occurrence of preformed pairs above the superconducting transi- 
tion. This RVB concept explains the decrease of the uniform spin susceptibility 
and the reduction of the specific heat below T*. The holes (vacancies) are re- 
sponsible for transport in the ab plane; the conductivity spectral weight in 
the ab plane is given by the hole concentration and is therefore unaffected by 
singlet formation. For c-axis conductivity an electron is transported between 
planes; since an electron carries a spin-1/2, it is necessary to break a singlet. 
This qualitatively explains the strong reduction in c-axis transport at T*. 

This RVB concept provides a qualitative understanding of the pseudogap re- 
gion. At a more quantitative level, it has been shown that this RVB theory 
gives an excellent description of the key features of the superconducting state 
[11,12,13]. On the other hand, the situation for the normal phase (above T^) 
is less satisfactory; the RVB theory can be formulated in a gauge theory but 
then approximations are required to evaluate it [14]. In fact, the instability 
of the Landau-Fermi liquid appears to be more than just an instability to su- 
perconductivity. For example, unlike the superconducting energy gap, which 
vanishes at isolated points on the Fermi surface (the d-wave nodes), the pseu- 
dogap vanishes along segments of the Fermi surface centered at these nodes. 
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known as Fermi arcs [15,16,17,18]. This strange phase separation in momentum 
space is a key feature of the pseudogap phase. Nernst effects become apprecia- 
ble above T^, (but below a temperature T^get < T*) indicating the existence 
of preformed pairs above [19]. Additionally, as a result of the prominent 
Hubbard interaction U, Mott physics still plays a crucial role as observed in 
optical transport [20]. Therefore, the instability of the Landau- Fermi liquid in 
the pseudogap phase could be a precursor of the Mott insulating phase. 

This has been a source of motivation to search for (doped) Mott insulators of 
the spin- liquid sort. Generally, spin liquids are more common in low dimensions 
where quantum fluctuations can suppress magnetism. Quasi-one dimensional 
(ID) ladder materials [21,22] are promising in this regard and have received 
extensive attention, particularly the two-leg ladder [23]. Interacting fermions 
moving on a single chain form the celebrated Luttinger liquid. But when chains 
are assembled to form ladders this gives rise to remarkable changes in the 
physics. The two-leg (2-leg) ladder at half-filling is an example of insulating 
spin liquid that can be regarded as a form of short-range RVB state; there 
is a charge gap and a spin gap in the strong-coupling [24] as well as in the 
weak-coupling regime [25]. This insulating spin liquid phase can be viewed as 
a quantum disordered superconductor; it exhibits pairing, with an approxi- 
mate d-wave symmetry. This phase is then referred to as the D-Mott phase 
[26], which exhibits an enlarged S0(8) symmetry in the weak-coupling regime. 
Moreover, upon doping, the holes bind in relative d-wave pairs and the two-leg 
ladder exhibits quasi-long range superconducting (d-wave) pairing correlations 
[25,27]. This behavior is reminiscent of that seen in the underdoped cuprate 
superconductors. The glue for the pairing of holes in the two-leg ladder system 
may be viewed as the spin gap formation; the holes pair to avoid the breaking 
of the resonating singlets. On the other hand, it is instructive to notice that 
spinless fermions can also support unconventional (p-wave) superconductiv- 
ity in the weak-interaction regime [28]. This result shows that the interchain 
hopping not only favors single-particle hopping between chains, but it also 
strongly reinforces the Cooper channel in the ladder problem. 

The three-leg (3-leg) ladder is also very interesting. At half-filling, the strong 
coupling limit is equivalent to the 3-leg Heisenberg spin- 1/2 antiferromag- 
netic ladder. This has been well studied and reduces in the low-energy sector 
to an effective single-chain Heisenberg model with longer range (but unfrus- 
trated coupling) [29]. In fact, it is well established that even-leg ladders have 
a spin gap and odd-leg ladders have one gapless spinon mode [23]: this is 
the well-known odd-even effect for spin ladders. Interestingly, we have shown 
that the weak-coupling limit also shows a similar behavior [30]. Here, it is 
certainly worth mentioning that while for (undoped) spin ladders, theory and 
experiment are in agreement [31], superconductivity has up to now only been 
observed in a two- leg ladder material under high pressure [22]. Upon doping 
the 3-leg ladder, both in the strong [32,33] and weak-coupling [30] limit, one 
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finds that the holes enter the channel with odd parity under reflection about 
the central leg and form a Luttinger liquid, whereas the even parity channels 
remain at the stoichiometric filling and continue to form an insulating spin 
liquid (RVB state). The result is a finite region of hole doping where the orig- 
inal Fermi surface with 3 bands (or 6 Fermi points) is partially truncated to 2 
Fermi points, but without a broken translational symmetry. The key is Umk- 
lapp scattering which introduces a charge gap in even parity channels even 
away from half-filling. We have also shown that this truncation of the Fermi 
surface is robust to the inclusion of a small repulsive next-nearest neighbor 
hopping t' [34] which should be added to any realistic (1-band) description of 
the cuprate systems. The 3-leg ladder is a clear example of a partially trun- 
cated Fermi surface through the formation of an RVB state over part of the 
Fermi surface. However, in the slightly doped 3-leg ladder, the Fermi surface 
reduces to nodal points and not to Fermi arcs as in high-Tc cuprates. 

The ladders exhibit an odd-even effect whereas the 2D system has two gapless 
magnon modes. It is then interesting to know how the quasi-lD system evolves 
towards the 2D counterpart. In antiferromagnetic spin ladders, using the non- 
linear sigma model [35], one can show that the spin gap for even-leg ladders 
disappear as exp — A^, whereas at weak [/, the spin gap disappears double- 
exponentially as a function of the number of legs N [36]. The next question 
is: "what is the relation between antiferromagnetic fluctuations and super- 
conductivity?". The quasi-lD approach, in the limit of large number of legs 
A'^ — > -|-oo, also allows to rigorously show that short-range antiferromagnetic 
fluctuations stabilize d-wave superconductivity for not too large doping levels 
[36]. A Kohn-Luttinger type attraction [37] is generated by antiferromagnetic 
processes like in two 2D Hubbard models [38,39,40]. In two dimensions, renor- 
malization group (RG) methods have been developed to treat weak to strong 
interactions. The antiferromagnetic and d-wave superconducting phase can be 
predicted from the examination of the RG flow, e.g., with a Fermi surface con- 
taining the saddle points (tt, 0) and (0, vr) (and Fermi surface nesting as well as 
van Hove singularities) [41,42,43]. The theory has been pushed further through 
the functional renormalization groTip by Honerkamp et al. [44]. A spin-density 
wave analysis in the ordered state and a discussion on the evolution of the 
Fermi surface in underdoped cuprates can be found in Ref. [47]. 

Even though a proper (rigorous) theory of the pseudogap phase is still miss- 
ing, some progress has also been accomplished recently. First, the quasi- ID 
approach shows that one can still have a truncation of the Fermi surface when 
N — > -|-oo; the antinodal points form an insulating spin liquid similar to the 2- 
leg ladder whereas the nodal regions are gapless. Moreover, this allows to prove 
that the antiferromagnetic (RVB) fluctuations (in the antinodal directions) 
will reinforce d-wave superconductivity at low energy in the vicinity of the 
nodal directions, resulting in d-wave superconductivity in the nodal regions. 
Second, by including a small repulsive next-nearest neighbor hopping t' , RG 
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arguments show that a similar scenario can emerge in two dimensions close to 
half-filling due to the effect of umklapp scattering [45]; this requires an interac- 
tion strength larger than Uc where Uc/t oc ln~^(t/t'). Those approaches suggest 
the appearance of two energy scales close to the Mott state: T* embodies the 
energy scale at which umklapp scatterings still play a crucial role in the antin- 
odal regions due to the proximity of the Mott state, producing an (insulating) 
RVB liquid, and Tc represents the superconducting energy scale at which the 
Fermi arcs become superconducting as a result of the prominent antiferro- 
magnetic fluctuations in the antinodal regions. However, a proper low-energy 
theory is lacking and has only been explored through numerical techniques 
[46]. The main difficulty is to build an unbiased 2D low-energy theory that 
include all the (divergent) particle-particle and particle-hole channels, namely, 
forward, backward (Cooper), umklapp, and antiferromagnetic processes. On 
the other hand, a phenomenological theory of the pseudogap phase [48,49] can 
be built by analogy with the ladder systems and an array of coupled ladders 
[50]. The proposal of a Fermi surface consisting of 4 disconnected arcs has 
parallels to gauge theory calculations for the doped t — J model [14,51]. 

Finally, it is also worth mentioning the recent development of Dynamical Mean 
Field Theory (DMFT) [52,53] and cluster type approaches [54]. 

Ultracold atoms in optical lattices are promising simulators of complex many- 
body Hamiltonians that arise in condensed-matter systems. They embody very 
clean systems which can be tuned in a very controlled manner from weak to 
strong coupling [55] . The Hubbard model has been recently realized with both 
repulsive and attractive interactions [56,57,58]. On the other hand, it has been 
shown for some time that the minimal 2D unit (the 4-site plaquette on the 
square lattice) can sustain d-wavc pairing physics [59]. The single plaquette 
has been the starting point of various approaches to the Hubbard model in 
two dimensions [60,61]. The hole pair that is created has a d-wave charac- 
ter and leads to d-wave superfluidity once the plaquettes are coupled. Cold 
atomic gases in optical lattices may allow to realize those "plaquette" systems 
[62] and then to prove that d-wave superfluidity can resist to inhomogcneous 
hopping amplitudes (even though one needs to improve the cooling procedure 
of fermions to access the Neel as well as the d-wave superconducting state) . 

The outline of this review article is as follows. 

First, we focus on ladder systems that have a continuous crossover between 
weak and strong interactions and show many salient features such as RVB 
physics and truncated Fermi surface as well as d-wave superconductivity close 
to the Mott state for spinful electrons. For weak interactions, we underline 
the emergence of an enlarged symmetry in the low energy sector. Then, we 
discuss the dimensional crossovers to two dimensions and discuss analytical 
expressions for the antiferromagnetic Hamiltonian as well as the d-wave su- 
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perconductor appearing upon doping. For large enough doping, one recovers 
a Fermi liquid when antiferromagnetic correlations disappear completely. We 
also describe the pseudogap phase emerging from the quasi-lD approach and 
show that at low energy the nodal regions are unstable towards d-wave super- 
conductivity as a result of antiferromagnetic fluctuations (or the emergence of 
RVB physics in the antinodal directions). We continue by reviewing the RG 
results in two dimensions and mean-field approaches for the strong-coupling 
doped t — J model. In particular, introducing a Gutzwiller projected BCS 
wavefunction and resorting to mean-field theory to solve the t — J model, 
gives a relatively good description of the d-wave superconducting phase. We 
also suggest that the pseudogap phase might be understood by drawing a par- 
allel with the A^-leg Hubbard ladder in the weak-coupling regime. We sum- 
marize the theoretical results and draw a comparison with experiment, both 
in the superconducting phase and in the pseudogap phase. We also address 
the crossover to the overdoped limit and the origin of the breakdown of the 
Landau-Fermi liquid near the onset of superconductivity. On the other hand, 
we show that cold atomic systems in optical lattices allow to probe several 
features of Hubbard models. In particular, we investigate the d-wave pairing 
in 2D plaquette Hamiltonians which can be implemented in optical lattices. 

In all these Hubbard models — from homogeneous (Hubbard model in two 
dimensions) to inhomogeneous models (plaquette and ladder models) — an- 
tiferromagnetic fluctuations favor a d-wave superconducting state. 

In fact, the occurrence of d-wave superconductivity in electron systems de- 
veloping prominent spin fluctuations at the wavevector Q = (vr, vr) (in two 
dimensions) can be intuited as follows. First consider the BCS gap equation, 

^scip) = -T.Vs{p-p')^^, (2) 
p' p' 

where Vs embodies the singlet channel interactiorQ (which is modified by 
antiferromagnetic spin fluctuations) [63]. When the interaction is attractive 
{Vs < 0) and slowly varying over the Fermi surface, one can see that this 
inevitably results in s-wave superconductivity, i.e., the gap has the same sign 
over all the Fermi surface. Now, let us consider repulsive interactions. Spin 
fluctuations will make the singlet channel interaction more positive at large 
momentum transfer p—p' = Q = (vr, vr). Now, suppose Asc{p') is positive near 
(0, vr) then this produces a negative gap at (7r,0), so that a dx2-y2 gap devel- 
ops. As shown below, this comes about from band-structure nesting effects in 
two dimensions close to half-filling or because of strong-coupling, short-range 
valence bond correlations such as in the (2-leg) Hubbard or t — J ladders. It 

/S.SC represents the superconducting gap and = (A|(-^(p) + ^^)^/^ (^p- embodies 
the kinetic term where we have subtracted the chemical potential). 
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is also interesting to note that in 1986 (approximately at the same time than 
the discovery of high-Tc cupratcs), Emery suggested that back-scattering from 
spin fluctuations could lead to the pairing of holes in the Bechgaard salts [64] . 
In addition, three papers in 1986 argued that antiferromagnetic spin fluctua- 
tions might mediate d-wave pairing in heavy fermion materials [38,39,65]. 



2 The Ladder prototype 

Even though the Hubbard chain does not support unconventional supercon- 
ductivity for repulsive interactions, the two-leg ladder supports unconventional 
superconductivity (p-wave superconductivity for spinless electrons [28] and d- 
wave superconductivity for spinful electrons). First, we give a justification to 
the occurrence of RVB physics and d-wave superconductivity in the two-leg 
spin (or t — J) ladder. Then, we discuss the weak-coupling limit and introduce 
the D-Mott phase, that is very similar to the insulating RVB state emerging 
for large interactions, but it exhibits a large SO (8) symmetry [26]. 

2. 1 From RVB physics to exotic superconductivity 

A pedagogical way to understand the emergence of RVB physics and d-wave 
superconductivity in ladders is through the large onsite interaction limit (or in 
the weak-interaction limit when the Mott gap in the chains exceeds the inter- 
chain tunneling amplitude [66]). At half-filling, the system then is equivalent 
to a two-leg spin ladder described by a superexchange J = 4t'^/U > along 
the chains and an intrachain superexchange J± = At^/U > [23]. Let us start 
with the limit J± ^ J. To minimize energy, the spins of the two chains will 
form singlet bonds, resulting in an RVB-type wavefunction: 

\RVB > = Pn=oEPi'-j)i4iAjl - dlAj^lVac >, (3) 

where di(2)is annihilates an electron in chain i of the ladder, \Vac) refers to 
the vacuum, and Pd=o is the Gutzwiller's projector [67] leaving only singly- 
occupied sites with spins, 

Pd=o = I[ n (l-<T<i)- (4) 

i a=l,2 

By construction, the function F{i—j) decays exponentially for distances larger 
than ^ cx 1/J±. From this simple analysis, one can conclude that in the limit 
J_L > J, the two-leg ladder has a Mott gap of the order of U and a spin gap 
~ J±. Additionally, the (gapped) S — 1 magnon excitation, corresponding to 
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Fig. 2. Two-leg ladder: the holes pair to minimize the cost of magnetic energy. 

break a singlet, has a dispersion of the form [68] (of course, the magnon is 
massive due to the prominent spin gap ~ J_l): 

uj{k) = J± + Jcos{k) + —-(3 - cos(2A;)) + ... . (5) 

4J_L 

Here, k refers to the wavevector associated with the direction along the chains, 
and for simplicity, the lattice spacing will be fixed to 1. The minimum of energy 
is at A; = vr emphasizing the antiferromagnetic correlations along the chains. 
To see the connection with d-wave superconductivity, it is useful to rewrite this 
rung-singlet state in terms of the bonding (1) and antibonding (2) operators: 

dl/2is = ^ (^li. ± ^2i.) . (6) 

Since the function F decays exponentially with space, one may approximate 
\RVB) by the state (the electron spins across the rungs of the ladder are 
locked into singlets): 

\RVB > ~ E(*LT*Li - "^lAWac > . (7) 

i 

This is equivalent to add a singlet (Cooper) pair into the bonding and anti- 
bonding orbitals. This paired form is suggestive of d-wave superconductivity. 
When viewed in momentum space, the ground state of a superconductor is 
a product of singlet pairs with zero center of mass momentum at different 
points around the Fermi surface. In an s-wave superconductor the pairs are 
all added with the same sign, but if the pairs are formed with a relative angu- 
lar momentum {e.g., d-wave) sign changes are expected. The connection with 
d-wave superconductivity will be made more concrete in the next subsection 
when representing the bonding and antibonding bands in the two-dimensional 
Brillouin zone. To conclude this analysis at large J±/J ^ 1, when equally 
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doping the two chains of the laddei[fj, one expects the emergence of d-wave 
superconductivity in the system to minimize the magnetic energy (see Fig. 2). 

It is also instructive to observe that close to half-filling superconducting cor- 
relations will exhibit a universal form. In the limit of large J± one may define 
the superconducting order parameter as: 

^^ = {d\,A^l-d\ildU)■ (8) 

Owing to the fact that in one dimension, one cannot have a true symmetry 
breaking, in the d-wave superconducting phase one expects: 

(AjA,)oc|z-jr'^", (9) 

where r]scd ^ 2 (when superconducting fluctuations are prominent). Close to 
half-filling, the exponent rjscd can be intuited as follows [70]. The low-energy 
Hilbert space consists of rungs of the ladder where either both sites are occu- 
pied and form a singlet or both sites are empty. Singly occupied sites or rungs 
with a triplet lie higher by an energy of order • It is convenient to define the 
state where all rungs are occupied as the vaccum such that an empty rung can 
be visualized as the creation of a (hard-core) boson h\ = Aj. In second order 
perturbation theory in t and J then one obtains the low-energy Hamiltonian: 

Heff = -teff Y^{h\hi+i + h.C.) + VeffY^ riiUi+i, (10) 
i i 

where = h\hi < 1, tg// = 8tV(3Jx), and K// = (16tV3 - 3JV8)/ J^. In the 
dilute limit for the hole pairs (close to half-filling) one can neglect the effect of 
the intersite interaction Veff and the Hamiltonian turns into a model of free 
fermions through the Jordan- Wigner transformation. By analogy to the XY 
chain, then one predicts the universal correlations [70]: 

(AlA,)cxK-jr^/'- (11) 

This power law close to half-filling (and the RVB picture at half-filling) appears 
for all values of J±/ J ^ [71]. (Remember that superconductivity has up to 
now only been observed in a two-leg ladder material under high pressure [22].) 

2.2 Weak- coupling regime 

Below, we demonstrate thoroughly that the physics still remains the same 
in the weak-coupling regime U <^ (t±,t). We apply KG techniques in the 

^ Here, we consider homogeneous doping in the two chains (below, 5 represents the 
hole density per site). On the other hand, a substantial asymmetry in the chemical 
potentials of the two chains is well-known to be pair-breaking [69] . 
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Fig. 3. Band structure of the two-leg ladder. 

interaction U along the lines of our Ref. [28]. In this limit, it is a good approach 
to diagonalize the kinetic term first; going to the momentum space, one gets 
two bands (the bonding and antibonding bands; see Eq. (6)). The kinetic term 
takes the precise form: 

Ho=Y.f dkE,{k)^l{k)^,{k), (12) 

4 = 1,2"^ 

where the dispersion relations obey: 

Ei{k) = Tt±-2tcos{k). (13) 

Band 1 is the bonding band and band 2 the antibonding band. See fig. 3. 

By analogy with the two-dimensional case, the associated transverse momenta 
are denoted k± = (0, vrji] (consult Fig. 4). Since we focus on the low-energy 
physics, we linearize the dispersion relation Ei around the Fermi momenta 
^kpi, which are fixed by the chemical potential /i = Ei{kFi) and the filling n 
corresponding to the electron density per site. At half-filling Ei{kFi) = 0. 

We first study the spinless fermion case where {kpi + kF2) = 2n7r and the RG 
equations can be solved analytically [28]. The Fermi wavevectors are given by: 

kpi = vrn ± arcsin f -, — - | . (14) 

Utsin vm / ^ ' 



The Fermi velocities thus obey: 



Vi-V2 = - r. 15 

tan(7rn) 



^ For the 2-leg ladder, we choose the convention E-i (k) = —t± cos k± — 2t cos k where 
= for the bonding band and k± = n for the antibonding band. This allows us 
to make an analogy with the antinodal regions of the 2D half-filled Fermi surface. 
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It should be noted that for spinlcss ladders, the filhng is < n < 1 and 
the hole doping (hole density per site) away from half-filling corresponds to 
5 = 0.5 — 11. Now, the effect of the interchain hopping is totally included in 
the two velocities Vi and ^2 which are only equivalent at half-filling. In this 
subsection, we do not consider the half-filled case, {kpi + kpT) — tt and focus 
on the system away from the insulating regime where the difference in the 
velocities (or a finite t±_ value) has the remarkable effect of driving the system 
to a superconducting state for repulsive interactions t/ > 0. 

Including all interactions allowed by symmetry (leaving away completely chiral 
interactions that do not affect the fixed point properties; see Appendix A), in 
momentum space, the Hamiltonian takes the form H — Hq-\- Hint where 

Ho^^Vi [ dkk (^Uk)^,Rik) - ^L(fc)^a(fc)) , (16) 

1=1,2 

and 



Hint — J dkidk'2.dk^dki5{ki + ks — k2 — k^) (17) 

X [ci^l^(/Ci)^iR(fc2)*Ii(fc3)*lL(fc4) +C2(1 ^ 2) 
+ /l2(^lK(fci)^li^(fc2)^L(^3)^2L(M + 1^2) 

+ ci2(^l«(fci)*2i?(A;2)*L(A;3)*2L(A;4) + 1 ^ 2)]. 

Here, /12 describes the interband forward scattering, C12 the interband Cooper 
coupling, and q the intraband Cooper coupling. Keep in mind that in this 
subsection we drop out all the umklapp terms (see Fig. 4) because we assume 
that the system is away from half-filling. 

For spinless fermions, the bare values are given by (Appendix A): 

ci = C2 = /12 = C12 = [/ > 0. (18) 

Those Cooper and forward interaction channels are reminiscent of the coupling 
channels in a 2-patch model in two dimensions; see Fig. 4 and Fig. 10. 

To second order in the interaction [/, the perturbative RG equations read [28]: 
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Fig. 4. Channels of interaction for the two-leg ladder in momentum space. A two-leg 
ladder produces 4 Fermi points in momentum space (it should be noticed that in 
this figure we have chosen kpi and kF2 arbitrarily but we have set k± = for 
the bonding band and k± = ir for the antibonding band.) In this basis, d-wave 
superconductivity emerges because the superconducting order parameter changes 
of sign when going from one band to the other (or under a 7r/2 rotation). 

dC2 1 2 

dl 2'KVi 
dh2 ^ 1 2 

dl ~ Tt{Vi + V2)^^'^ 
dCi2_Ci2 ( 2/i2 Cl 
dl TC \Vi + V2 2vi 

Those RG equations can be found quite easily using the field theory rules 
of Eq. (A. 13). The energy scale (temperature scale) is related to / through 
E ~ te~'. The plus and minus signs result from particle- hole respectively and 
particle-particle diagrams. For specific details, consult Appendix A. 



Keeping only the terms quadratic in the coupling constants, the set of 4 dif- 
ferential equations turn into one differential equation for f^. First, using all 
the RG equations, it is relatively easy to derive: 



dci2 

~dr 



Cl2 



12 



{Vl+V2)U' 
7rViV2 



(20) 
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where 

+ >0. (21) 

27rt'iW2(wi + V2) 

Then, resorting to the third equation, one can show that, 

j^-{h2-BUf + CU', (22) 

where 



B: 



2(vi + V2) 



> 



—vf + 6vfv2 + 6f I + 6V1V2 — V2 

{AViV2 + {V^ + V2yf 



(23) 



The solution of Eq. (22) is different for C < and C > 0. For C > 0, all 
the couplings diverge whereas for C < 0, C12 flows to zero and all the other 
couplings remain small (~ U). In terms of the velocities, C > means: 

1/7 < V1/V2 < 7. (24) 

One can check that the weak-coupling fixed point appearing for C < is 
stable at all the orders because whatever the order is, the RG equations are 
always multiplied at least once with C12. 

The low-energy physics can be analyzed thoroughly through standard bosoniza- 
tion techniques [71,72]. The Hamiltonian takes the form [28]: 



^-E/<^(f + ^)(9.«^ + (|-^)n? (26) 

+ Jdx^ {dAidx(t>2 - nin2) - cos (V4^(^i - ^2)) . 

(Again, the lattice spacing is fixed to 1 and the Klein factors ensuring the cor- 
rect anticommutation rules for the fermions have been chosen as r]iRrjiLr)2R'r)2L — 
1.) Here, dx4>i represents the charge fluctuations in each band and 6i corre- 
sponds to the conjugate superfluid phase; we have introduced Ilj = dxOi. A 
flow to strong coupling of C12 (quasiclassically) results in the "pinning" of the 
superfluid phase ^1 — ^2 = in order to minimize the energy, and in a single 
gapless mode, while for C12 — > 0, two gapless modes are present. 

A superconducting phase is in fact stabilized when C12 flows to strong cou- 
plings, i.e., for Vi/v2 < 7. The key point is to resort to the canonical transfor- 
mations (f)± = (01 ± 02) /\/2 and n± = (Hi ± n2) / v^. The mode ^_ is pinned 
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to zero due to the prominent term Cu whereas the symmetric mode 0+ and 
6^ remain gapless. The effective Hamiltonian at low energy takes the form, 

H+ = Jdx^ 1^ J_(40^)2 + K+ul^ + A9,0+9,.0_, (26) 

We have ignored a term of the form 11+ n_ because this term is irrelevant 
when Ci2 flows to strong couplings. We identify: 



v± 




(27) 



27r(t;i + V2) - (ci + C2± 2f: 
^ 2tt{vi + V2) + (ci + C2 ± 2/] 



A: 



2 47r ■ 



We have introduced the velocity and the Luttinger parameter of the antisym- 
metric bosonic mode because they will be used belo43. 



Now, we discuss the main correlation functions in the problem. For repulsive 
interactions, the charge density and superconducting pairing fluctuations with 
the most divergent susceptibilities are 



5. 



ocexp (i\/2TT(l)^^ cos (^VQjcOJ^ 



(2J 



and 



Osc = diRd2L + d2RdiL = ^IR^lL " ^2R^2L (29) 

ocexp (j.V2tt6^^ cos (^V27t6_^ , 

where dn/R represent the annihilation operators for the fermions in chain 
i. Note that interchain superconductivity means intraband superconductivity. 
Additionally, for spinless fermions, the superconducting order parameter must 
be odd under parity, Osc{—x) = —Osc{x), which is reminiscent of p-wave like 
superconductivity. The cdw and sc correlation functions are given by: 

^ At the fixed point, as a result of /12 > and q < 0, A'_ > that will favor the 
pinning of the superfluid phase difference Oi — 62. 

^ At half-filling, for spinless fermions, the ground state is a charge density wave 
with density alternating between the chains (n = 1/2). The symmetric phase 0+ 
becomes disordered due to the umklapp scatterings. 
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Fig. 5. Our phase diagram of the spinless ladder for tj^/t = 1 and U/t = 0.2 [28]. 
Tq ^ t represents a high-energy cutoff. The solid line shows the crossover to the 
(p-wave) superconducting phase where phase coherence between the bands arises 
(close to half-filling, this occurs at the temperature scale Tc ~ te~'^'^^/^) and the 
dashed line indicates when superconducting fluctuations appear in band 1. 

{OlaSx)0^^m oc x--^ (30) 
{Ol{x)Osom oc x-"\ 

where [28]: 

^ = 1 K^K.y^r (31) 

In the superconducting phase, to leading order in t^jt and in the doping 
(5 = 0.5 — n, one can approximate [28]: 

7-l4(^i)^*^. (32) 

It is important to observe that as a result of the finite interchain hopping 
oc A leading to 7 > 1, the system will develop superconducting pairing 
correlation functions when ^1/^2 < 7 (for = t, this means < 6 < 0.22). 

For Vi/v2 > 7, the fixed points of the couplings are such that C12 = 0, and Ci, 
C2, and /12 arc of the order of U. It is interesting to observe that in a short 
range of ratios, 7 < Vi/v2 < 8, the low temperature fixed point is such that 
Ci < (as a result of the small velocity V2) whereas C2 > 0, implying that the 
bonding band 1 will still develop prominent superconducting fluctuations while 
the (antibonding) band 2 will exhibit charge density wave fluctuations. In this 
case, the superconducting correlations occur both between the chains as well as 
along the chains. This phase can be seen as a precursor of the superconducting 
phase occurring at Vi/v2 < 7: the system develop phase fluctuations but no 
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phase coherence between the bands. Finally, for ratios fi/t'2 > 8. the couplings 
Cj > and as a result in both bands charge density wave fluctuations dominate. 
This corresponds to the usual metallic phase in ID: the Luttinger phase. 

Our phase diagram of the spinless ladder is shown in Fig. 5 [28] . 

This analysis can also be useful for other two-band systems such as carbon 
nanotubes on a superconducting substrate [73] or quantum wires [74]. Now, 
we extend the analysis to the case of fermions with spin following Ref. [27], 
and show that the system will develop superconducting d-wave fluctuations 
as well as a prominent spin gap, as a reminiscence of the large U analysis. 

Away from half-fllling the interaction channels in the band basis take the 
compact form: 



Hint — j dx {jijJmiJLjj - ftjJRiiJhjj + cljjRijJuj - c^ijJmjJLij) (33) 
+ dx (c'^iJRiiJLii — clj^JmiJui) ■ 

1=1,2'' 

The only diffference with the spinless case is that now the system develops 
Cooper and forward scattering channels in the spin sector as well. For a pre- 
cise relation between the c^j and parameters and the usual g-ology of the 
Luttinger theory, consult Appendix A. We have exploited the U(l) symme- 
try for charge and SU(2) symmetry for spin, introducing the charge and spin 
densities: 

Jpij = E ^U^.Ps and jr^ = Iy^ ^U^'^Pi.', (34) 

s ^ s,s' 

where p — {L,R) denotes the direction (chirality), i represents the band, 
denote the Pauli matrices, and s the spin. By symmetry, one gets c^j* = Cjf 
and fif — fjf. Here, the Hubbard bare values are given by (see Appendix A): 

4i- = = 44 = 4 = (35) 



In contrast to the spinless case, the (one-loop) RG equations cannot be solved 
analytically, such that one needs to perform a numerical integration. On the 
other hand, a characteristic property of the N-leg ladders is that the couplings 
flow towards universal ratios as the initial value U/t is reduced; in particular, 
away from half-flUing and at energies smaller than E ~ te~*^^, one gets [27]: 

t ~ 4c^2 = c'i2, cljvi = - -1, and fl^ ~ 0. (36) 

One can bosonize the Hamiltonian; the low-energy flxed point is governed by: 
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Hint = J dx - \cl^ \ cos{V8TT(f)is) - jcgal cos{V8Tr(f)2s) (37) 
- 20^2 cos{VAtt6c^) (cos{V^(f)s-) + cos(V47r0s+)) • 

Minimizing the energy leads to the pinning of the spin modes in each band 
(pis = (p2s = (or (j)s+ = {(pis + (t>2s)/V^ = and = {(pu - 02s)/\/2 = 0) 
and of the antisymmetric charge mode 9c- similar to the spinless case ^ I. Since 
both spin modes are pinned, this shows that the ladder has a spin gap when 
U is small. 

The total charge mode remains gapless similar to the spinless case and one 
also finds that the superconducting pairing correlation function is the most 
dominant one. The (singlet) pair-field operator reads: 

A, = ^,7?T^iLi + *iLT*ii?i « VnVjie''^''^ cos( v^0,-,) ; (38) 

we have introduced the Klein factors rjjg (ensuring the correct anticommuta- 
tion rules between two fermionic operators) [71,72] such that ^7it'7ii'72T^2i = 1- 

One obtains that the pairing correlation function decays as: 

(A|(x)Aj(0)) oc x~ 

i/{2i^c+)^ (39) 

where Kc+ ~ 1 is the Luttinger parameter of the total charge sector. We check 
that the pairing correlations in the spin-1/2 case decay as oc x"^/^, whereas 
all the other correlation functions decay approximately as The spin-1/2 
two-leg ladder is superconducting both for t± > U and for t± < U, whereas 
for the spinless ladder superconductivity disappears when ^ or A ^ 0. 

Furthermore, for spinful electrons, one can check that the superconductiv- 
ity has indeed an approximate d-wave symmetry meaning that the pair-field 
operator has a different sign in band 1 and in band 2, i.e., (AIA2) < 0. 



2. 3 D-Mott state from weak coupling 

Similar to the t — J ladder, the two-leg Hubbard ladder at half-filling is an 
insulating spin liquid or a Mott insulator with a spin gap. In fact, the D-Mott 
insulator [26] is a concrete example of disordered d-wave superconductor: the 
phase coherence between the bands is already present, but the total charge 
mode is gapped due to the presence of extra umklapp scattering processes. 



" In particular, the fixed point leads to the following Luttinger parameters: Ks_ < 
1, Ks+ > 0, and Kc- > 1. 
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Fig. 6. Coupling ratios in the half- filled two- leg Hubbard ladder versus the next 
nearest-neighbor hopping —t' following the analysis of our Ref. [34]. For a reasonable 
value of t' ~ 0.1 — 0.2t, the system remains insulating and still sustains a spin gap. 

More precisely, since kpi + kp2 = vr at half- filling, this allows for the additional 
two-particle umklapp processes (for an illustration, see Fig. 4): 



H 



Int,hf 



= Jdx (^u522l(4l2^L21 + h.C.) - U^i22l((4l2)^ " /l21 + h.C.) (40) 
+ <122(4ll^L22 + 422^^11 + h.C.)y 



We have introduced the umklapp operators: 



(41) 



where e = —ia^ and Ipij = Ipji, and: 



^pij 2 P^^^ Jss'^pjs 

s,s' 



(42) 



The initial bare values are given by (consult Appendix A, i.e., Eq. (A. 19)): 



'^'^kkkk = ^"^IkVk = U and -uf rr, = 0. 



(43) 



The D-Mott phase is dominated by the following set of couplings (see Fig. 6): 



t ~ 4c^2 = 4/{=2 = cl, = -4 



4m? 



1221 



'^1122 ~ ^1221' 



(44) 



and = ~ /12 ~ 0. Both the Cooper and elastic umklapp channels flow 
to strong couplings at half-filling; therefore, to properly describe the low-energy 
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fixed point, one needs an exact strong-coupling theory that can distinguish be- 
tween superconductivity or insulating phase. In one dimension, one can apply 
bosonization [71,72]. The bosonized form of the umklapp interactions is: 



Since ^5221 — '*^i22i/4 in the D-Mott phase, one verifies that the system has 
still a spin gap: both 0^- and 0s+ are pinned to zero. Additionally, the phase 
coherence between bands is already present through the pinning of the anti- 
symmetric superfluid mode 6^-. On the other hand, at half-filling, the total 
charge mode (j)c+ is also pinned that makes the system insulating: when the 
mode 9c- is pinned, the current density only involves the total charge mode 
j — Vc+Kc+Ilc+, and half-filling when 0c+ is pinned this implies j — 0. In the 
weak-coupling limit, the charge gap and the spin gap are equal (~ te~'"^/^). 

The pair field phase 9c+ being the conjugate field will fluctuate wildly. These 
quantum fluctuations will destroy the power-law ID superconducting phase, 
leading to an exponential decaying pair-field correlation function. In fact, fur- 
ther progress can be accomplished by re-fermionizing the low-energy Hamil- 
tonian [26]: the Hamiltonian turns into an SO (8) Gross-Neveu model that has 
been well studied by particle field theorists. The S0(8) Gross-Neveu has a 
peculiar symmetry (triality [75]) that allows to access the energies of various 
excited states. In particular, it has been shown that the energy of the lowest 
excited state with the quantum numbers of an electron is equal to the energy 
of the Cooper pair. This demonstrates pairing in the D-Mott phase. 

Finally, it is also important to underline that the charge gap and the spin 
gap are robust if one includes a moderately small (frustrated) next-nearest 
neighbor —t' hopping with t' > (relevant in high- Tc cuprates); see Fig. 6 and 
our Ref. [34]. In particular, for a relatively small t' ~ 0.1 — 0.2t, the RG flow 
indicates that the system still exhibits a large symmetry and only the couplings 
cli and C22 progressively deviate from the other (diverging) couplings. 

2.4 Emergent symmetries 

In Sec. 2.3, we have focused on properties of the D-Mott phase occurring in the 
weak-coupling regime of the two-leg Hubbard ladder; it exhibits a remarkable 
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S0(8) symmetry [26]. Interestingly, there exists an S0(5) subalgcbra of the 
S0(8) group whose vector representation unifies superconductivity and anti- 
ferromagnetism [26]. Thus the SO (5) symmetry, proposed by Zhang [76] as a 
phenomenological model for the cuprates, is shared by the D-Mott phase (and 
by the d-wave superconducting phase in the Hubbard ladder that maintains 
an S0(6) symmetry [77]). Below, we discuss the S0(5) theory. 

Let us consider two sites and define the fermion operators c and d on those 
two sites. The antiferromagnetic order parameter is defined as: 

= -(cV'^c-dVd), (46) 

a" represent the Pauli matrices acting on the spin space and the 'intersite' 
superconducting order parameter reads: 

At = -icW^ l(_,t^t + ,t^t). (47) 

One can group these five components and form a vector: (rii, ^2, ns, n4, ns) 
where rii = (A^ + A)/2, n2 = m^, = m^, n^^ = m^, and ns = (A''' — A)/2i 
called the superspin (it contains both superconducting and antiferromagnetic 
spin components). The S0(5) theory is inspired from the Landau- Ginzburg 
theory; on the other hand, the S0(5) theory extends the Landau-Ginzburg 
approach in several ways. First, it assumes an approximately symmetric inter- 
action potential between the antiferromagnetism and superconducting phases 
(in the underdoped regime of the cuprates). Second, it introduces a full set of 
dynamically variables canonically conjugate to the superspin (including the 
total spin, the total charge, and the so called tt operators that perform ro- 
tations from antiferromagnetism to superconductivity and vice versa) [78]. 
Through the introduction of the dynamically conjugate variables, the SO (5) 
theory is capable of describing quantum disordered phases. 

Microscopic SO (5) models can be explicitly constructed from ladder or bilayer 
systems [79]. Introducing the two sites c and d (for the ladder: c belongs to a 
rung and d to the other rung), the most general interaction Hamiltonian is: 



H2sites = U (net - ^) (u,^ - + {c ^ d) + V{n, - 1) {ud - 1) (48) 
+ JSdSc-i^{nc + nd), 

where ric = J2a ''^ca ^-nd a = t or J,. In fact, an important condition ensures the 
local S0(5) symmetry: J = 4{U + V). Lin, Balents, and Fisher have explored 
the possibility of SO (5) low-energy fixed points in the weak-coupling regime 
of the two- leg Hubbard ladder [26] ; when interactions are repulsive they found 
that the SO (5) system falls into the basin of attraction of the D-Mott phase. 
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On the other hand, the behavior of the spin-spectral function is different in 
the ladder system and in the SO (5) theory. The SO (5) theory predicts a sharp 
resonance at energy 2fi and momentum (vr, vr), the so-called vr resonance (orig- 
inally introduced by Zhang [76] to explain the 42meV neutron scattering peak 
in the superconducting cuprates). Lin, Balents, and Fisher have shown that, 
in addition to SO (5) symmetry, the delta-function vr resonance requires a non- 
zero condensate density in the superconducting phase. Since condensation is 
not possible in one dimension, this precludes a delta-function vr-resonance 
(however a weaker algebraic singularity is not precluded [26]). 

Regardless of the behavior in the vicinity of u; = 2/i, one expects spectral 
weight at energies below 2/i but above the spin gap. 

Interestingly, the optical conductivity at hal-filling exhibits a high-frequency 
preformed pair continuum as well as a sharp excitonic peak below the gap as 
a signature of the enlarged S0(8) symmetry [26,80]. 



3 Quasi-ID solution and Truncation of the Fermi surface 



Interestingly, the weak-coupling RG approach close and at half-filling can be 
extended to N bands with N > 2 [30]. The relation between the chain i and 
the band j reads: 

d^s{x) = E /^sin (^) ^,^{x). (49) 

We have exploited the open boundary conditions in the direction perpendicular 
to the chains (legs), such that the transverse eigenfunctions are standing waves. 
This brings the Hamiltonian into a diagonal form in momentum space, and 
for a number of legs N > 2, the dispersion relation of a given band j is: 

Ej{k) = -2tcosk- 2t^cos{nj/{N +1)) . (50) 



The Fermi momenta in each band kpj are determined by the chemical potential 
fij = Ejikpj) and the filling n = 2{TrN)-^ J2j kpj- At half-filling (yU = 0): 

kpi = vr — arccos ( — cos ( — . (51) 

^ \t \N+lJJ ^ ^ 

From the band structure, one can attribute a fictitious transverse momentum 
to each banclZ] ttj'/ (A^-|-1). It is relevant to note that for odd, there is always 



In fact, the transverse momenta are zizTTj/{N + 1). 
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a band, namely j = {N + l)/2 that lies in the nodal direction at (vr/2, vr/2) on 
the 2D half-filled Fermi surface whereas for N even all bands lie away from the 
nodal points. For an illustration, see Fig. 7. One can check that at half-filling 
and for t_i_ = t, all the Fermi points belong to the 2D umklapp surface. 

Again, since we are only interested in the low-energy physics, we linearize the 
dispersion relation at the Fermi surface, resulting in the Fermi velocities Vj = 
2tsm{kpj). At half-filling, one gets the important equalities (the lowerscript k 
referring to a given band k should not be confused with the wave- vector k) : 

Vf, = v-k = 2^Jt^-tl cos(7rA;/(iV+ 1))2, (52) 

where: 

k = N + l-k. (53) 
In the following, we will denote the bands {k, k) band pairs. 

It is important to note that the Fermi velocities fulfill: 

vi = vn < V2 = vn_i < ... . (54) 

The velocities are always maximum for bands in the vicinity of the nodal 
directions. For N not too large, these values of the velocities will lead to a 
hierarchy of energy scales, as shown below. Note also that at half-filling, the 
Fermi momenta of band pairs add up {kp^ + ^fs) = allowing interband 
umklapp processes to take place, similar to the two-leg Hubbard ladder. 

At half-filling, the interacting part of the Hamiltonian consists of forward. 
Cooper, and umklapp processes within a band and between different bands 
similar to the 2-leg ladder, as well as extra (non-)umklapp processes between 
4 bands (or 3 bands for an odd number of chains). On the other hand, when 
the number of chains is quite small, say, = 3, 4, one can rigorously show 
that the 4-band couplings (or 3-band couplings) remain arbitrary small at low 
energy; consult Appendix B and details in our Ref. [30]. This is an important 
characteristics of the small A^ limit (on the other hand, the 4-band couplings 
will become important in the limit of large A^ close to and at half-filling). 

Integrating the RG equation^, we find that the Fermi velocities Vk lead to a 
hierarchy of energy scales. El ~ te~°'''^^^^ with ~ 1, where the band pairs 
(A;, k) become successively frozen out and form a D-Mott state. For example, if 
one starts with the 3-leg ladder where, say, t±/t ^ a/2 such that f i = fs <^ f2, 
one gets the following rigorous results [30]. Firstly, the renormalization of 
single-band and two-band interaction channels involving the bands 1 and 3 is 
exactly the same as for the two-leg ladder and one gets the emergence of a 
D-Mott state at the energy scale ~ te~"^'^^/^. Second, there is an umklapp 

^ The RG equations are given in our Ref. [30]. 
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Fig. 7. The bands of the 3-leg (3-chain) Hubbard model for t± = t on the 2D 
umklapp surface. Band 2 hes in the nodal direction at (7r/2,7r/2). At half-filling, 
the system is insulating and is analogous to the Heisenberg chain at low energy. By 
slightly doping, one gets a truncated Fermi surface: the holes enter first the band 
with the lowest charge gap (i.e., band 2) whereas the bands 1 and 3 still form an 
insulating spin liquid (D-Mott state). When increasing the doping further, one gets 
a ID d-wave superconductor; band 2 then embodies the nodal (Luttinger) liquid. 

term for the band 2 since kp2 = 7r/2, and this opens a charge gap (but no 
spin gap) in band 2 at the energy scale E2 ~ te~"'^'"^^^ <^ E^. At half-filling, 
this allows us to conclude that the system has a gapless spin mode in perfect 
accordance with the odd-even effect for the spin ladders [23]. 

When doping, the first hole will enter into band 2 and the chemical potential 
jumps from to E^. Band 2 forms a (Luttinger) liquid whereas bands 1 and 
3 still form a D-Mott state and remain insulating ^ I: this demonstrates the 
possibility of a truncated Fermi surface close to half-filling in Hubbard-type 
models, as a result of the band structure which favors the hierarchy of velocities 
Vi = vj\f < f 2 = utv-i < ••• and then a non- uniform Mott gap along the two- 
dimensional Fermi surface. When increasing the doping further, as soon as the 
chemical potential will reach E^ then the bands 1 and 3 will form a d-wave 
superconductor whereas the band 2 will develop a nodal liquid (or a Luttinger 
liquid in one dimension). This fixed point is stable until very low energy scales 
because a nodal point cannot develop an Andreev process (where two paired 
electrons from the antinodal band 1 or 3 hop onto band 2 and produce a 
superconducting proximity effect in band 2) due to the d-wave symmetry [82]. 
Finally, the effect of doping on the 3-band system is illustrated in Fig. 7. 

For the half-filled 4-leg ladder, the bands 1 and 4 form a D-Mott state with a 
gap El and the bands 2 and 3 also form a D-Mott state with a smaller gap E2. 
The system can be viewed as a D-Mott state with a non-uniform Mott and spin 

^ This is consistent with the Lieb-Schultz-Mattis theorem [81]; when applied to 
the 3-leg ladder, this theorem predicts the existence of gapless excitations at a 
wavevector 3tt{1 — 6) away from half- filling. When holes enter only the "odd parity" 
channel (band 2 electron operator is odd under parity with respect to reflection 
about the middle chain) gapless excitations will have a wavevector 7r(l — 36). These 
two wavevectors are equivalent modulo (27r) so there is no inconsistency. 
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Fig. 8. When slightly doping the 3-leg ladder, the holes are situated on the outer 
legs (that means band 2). For the 4-leg ladder, the highest probability to find a spin 
singlet is when the two holes belong to the outer legs; the lowest probability to find 
a singlet is on the (inner) legs 2 and 3; consult Eq. (57) and our Ref. [30]. 

gap along the 2D half-filled Fermi surface. The whole system develops a spin 
gap at half-filling. As in 3-leg ladders, the channels decouple in pairs so that the 
energy to add a hole varies around the Fermi surface. Upon doping, the holes 
initially enter with the smallest energy gap. Therefore, the bands 2 and 3 first 
develop d-wave type superconducting correlations whereas the bands 1 and 
4 still develop a D-Mott gap. This is a concrete example where Mott physics 
can coexist with superconductivity. By increasing the doping, the system will 
develop a full d-wave superconducting gap and finally phase coherence between 
the band pairs will occur when adding a small Josephson coupling between 



band pain 



10 



Additionally, due to the phase coherence between band pairs i 



and j, one expects (here, we assume that interactions are very weak): 



(Aj(a;)A,(0))(xx- 



1/4 



(55) 



In fact, this successive opening of the Fermi surface with doping starting from 
the diagonals in the Brillouin zone will continue for a N-leg ladder system (for 
a discussion regarding the hmit N ^ +oo, see Sec. 4.2). 

Close to half-filling, it is also instructive to rewrite the band operators in terms 
of the chain operators. For the 3-leg ladder, one can check that the band 2 
annihilation electron operator is odd under parity (with respect to refiection 
about the middle chain) 

^2ps = ~y^{dips — dsps), (56) 



" For a large number of chains, the antiferromagnetic 4-band couplings automat- 
ically induce phase coherence between band pairs, i.e., d-wave superconductivity. 
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traducing that the unpaired holes will fill the outer legs of the ladder, as 
illustrated in Fig. 8. This is in accordance with the numerical results of Refs. 
[32,83] for the t — J Hamiltonian; it should be noted that in the case of strong 
interactions it is necessary to use numerical techniques to analyze the problem. 

For the 4-leg ladder close to half-filling, the superconducting order parameters 
of bands 2 and 3 are in accordance with d-wave pairing, (A2A3) ~ — 1, and 

A2 oc 0.22 {diR^d^n + d2B,^din + R ^ L) (57) 
+ 0.22 {dsR^d^Li + diR^dsLi + R ^ L) 

- 0.36 {diR^d^Li + d^Ri^diLi + R ^ L) 

- 0.14 {d2R^d^Li + d3R'fd2Li + R^ L). 

The singlets are on the top two legs, the bottom two legs, on the legs 1 and 
4, and with the lowest probability on the legs 2 and 3, similarly as found in 
Ref. [84] for the t — J model; consult Fig. 8. The 4-leg ladder in the strong 
interaction limit was also analyzed by Siller et al. [85]. The possibility of bi- 
pairing and stripe phase in the 4-leg ladder has been studied in detail [86]. 
Finally, the phases with N — 5,6 legs are in agreement with numerical works 
[84] , emphasizing the universal behavior emerging in the Hubbard ladders. 



4 Dimensional crossovers and Theories in two dimensions 



Now, we study thoroughly how the weakly interacting A^-leg Hubbard ladder 
evolves towards the two-dimensional — > +00. It is relevant to note 
that for t = t±, fi = fjv ~ 2Tit/N (consult Eq. (52)), leading to a singular 
behavior for large N (the quasi- ID analog of the van Hove singularities in 
two dimensions). A crucial difference between the quasi- ID and 2D case are 
the interactions which control the low-energy physics. Our quasi- ID approach 
is valid as long as the energy difference between two neighboring bands is 
larger than the largest energy scale in the system, i.e., te^*/^^. Using the band 
structure of Eq. (50) and considering t ~ t_L, this implies that: 

U ^t/\nN,tJ\nN. (58) 

Assuming this is fulfilled, one can take N ^ 00 and focus on the crossovers. 
Below, we thoroughly analyze the quasi-long range antiferromagnetic correla- 
tions emerging for the half-filled case and the d-wave superconducting state 
occurring away from half-filling. Close to half-filling, the situation is identical 
to the small limit: one gets a truncated Fermi surface and the antinodal 
directions form an (insulating) RVB state with preformed pairs. The RG equa- 
tions are a simple generalization of the 3-leg RG equations given in our Ref. 
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[30] and for the sake of clarity technical details are hidden in Appendix B. 
We will also discuss 2D approaches in the weak-coupling and strong-coupling 
regimes, which also support the RVB scenario close to half-filling (pseudogap 
phase). We like to emphasize that the d-wave superconducting state can be 
understood through a Gutzwiller-type projected BCS wave-function [12]. We 
discuss the electron's Green function in the pseudogap phase by analogy with 
the N leg ladder system [30] and with the array of coupled ladders [50]. We 
also address the crossover to the Fermi liquid regime taking place when the 
umkapp scattering (antiferromagnetic) processes vanish completely, making 
the Cooper channel irrelevant for purely repulsive interactions [87]. 

4-1 Antiferromagnetic long-range correlations 

For a finite and even number of legs (bands) A^, the ground state at half-filling 
is an insulating spin liquid characterized by the decoupling of band pairs {k, k). 
On the other hand, when increasing A^, since Vk — Vk+i ^ 1/N, this decoupling 
of band pairs will be suppressed and the "4-band" (antiferromagnetic) inter- 
actions of Fig. 9 will start to play a dominant role. Below, we show that this 
results in a uniform Mott gap and in long-range antiferromagnetic correlations. 

In fact, the forward spin coupling /^^ is prominent in the (quasi-lD) antifer- 
romagnetic phase. The sum of the couplings Rk = c^^ + /^^ then can be used 
to distinguish between the insulating D-Mott phase and the antiferromagnetic 
phase. More precisely, as shown earlier in the context of the two-leg ladder, 
in the insulating spin liquid (D-Mott) regime one expects that < as a 
result of the growing of c^^ (towards negative values) whereas in the antifer- 
romagnetic phase one finds that SU(2) single-band Cooper processes remain 
small (but they have become attractive) \cl.j^\ <^ t. Moreover, /^^ now grows 
in the antiferromagnetic phase explaining that becomes positive (see, for 
example, Ref. [36]). On the other hand, there is no broken symmetry in the 
quasi- ID antiferromagnetic phase and a spin liquid phase can still re-appear 
at half- filling at (very) low temperatures (below the antiferromagnetic phase). 

The energy scale at which the spin liq uid appears is precisely defined as the 
energy scale Ec{N) at which Ri ^ CQiJ. In fact, a rigorous analysis shows that 
the spin liquid only survives for very small energies smaller than [36]: 

Ec{N) ~ texp(-aexp(6A^)), (59) 

where a -C 1 and 6 is a constant of the order of unity (that depends on the 
ratio t±/t). Note that for the other band pairs {k, k), the corresponding energy 

The band pair (1, N) dominates in the RG process as the result of the van Hove 
singularity; vi = vn 2TTt/N. 
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Fig. 9. 4-band processes and D-wave Cooper channel in the large N limit. The square 
corresponds to the umklapp surface or the Fermi surface at half- filling {t = t±). 
Left: There are two types of antiferromagnetic (AFM) processes: umklapp processes 
(dashed arrows) and particle-hole processes (solid arrows). Right: The AFM pro- 
cesses which take place within a band pair are identical to the Cooper processes. 

scales are smaller since the velocities Vk are larger. For the spin ladder, the spin 
gap at large N decreases as, Jexp(— 0.68A^) [35]; while for large U the spin 
gap decreases exponentially, for small U the decrease is double- exponentially. 

In Appendix B.2, we study the phase for energies larger than Ec{N) and 
show that indeed it can be viewed as a precursor of the 2D antiferromagnetic 
phase. The antiferromagnetic type ground state occurs at an energy scale 

Eafm ^ te"^*/^, (60) 

where c is a function of t±/t; Eafm is enhanced for t± ^ t as a. result of the 
von Hove type singularity. Note that the 4-band antiferromagnetic processes 
of Fig. 9 with a weight oc are responsible for the RG instability at the 
energy scale Eafm- This also favors the occurrence of a uniform Mott gap. 
Moreover, the spin correlation functions now obey (for any A^, even or odd) 

{Si{x) ■ Si{0)) oc (-1)'^+' cos(7rx)/x^/^, (61) 

where the lowerscripts i and / represent the chains, and the space coordinate 
obeys x < 1/ Ec{N) since we are in the quasi-lD antiferromagnetic phase. 

4-2 Truncated Fermi surface and d-wave superconductivity 

By doping, the chemical potential /i couples to the total charge mode of each 
band pair (see Appendix B.2). Since all the 4-band interactions contain 

the field (pkic+i they will vanish. On the other hand, the effect of doping on the 
2-band interactions is the same as when doping the D-Mott state. Then, this 
will lead to a spin gap and d-wave like phase coherence between the bands k 
and k. The charge gaps close to the nodal point (7r/2, 7r/2) are the smallest ones 
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such that the holes enter there first. We observe that the effect of doping in this 
quasi-lD approach is very similar to decrease the number of chains (doping 
suppresses the antiferromagnetic 4-band processes at low energy scales). At 
very small doping, the bands close to the nodal directions will be slightly 
doped whereas the band pairs close to the antinodal directions will remain 
insulating as a result of the van Hove type singularities vi = vn ^ 2nt/N . It 
is instructive to observe that the limit of a large number of legs +00 
also supports the concept of a truncated Fermi surface close to half-filling. 

Now, we show that the nodal regions are in fact unstable towards d-wave 
superconductivity at low energy. At a very general level, one can prove that 
phase coherence between band pairs around the nodal directions will take place 
as soon as umklapp scattering and antiferromagnetic processes in those band 
pairs are cutoff by finite doping effect. Remember that close to half-fiUing the 
antinodal regions still form a D-Mott state because the charge gap is bigger 
as a result of the van Hove type singularity Vi — vn ^ 2Trt/N. 

Using the low-energy theory of Appendix B.2, at short distances, one can 
check that antiferromagnetism always coincides with superconductivity within 
a band pair (see Fig. 9, Right): 

- mfc • = -^'U^^LfeT^Lt^Kfei (+other terms), (62) 

and 

AlAfc = ^U^lLT^feLT^fefli (+other terms). (63) 
Now, by doping, one expects the superconducting fiuctuations within the band 
pair {k, k) in the nodal regions to become more prominent (by analogy to the 
2- leg ladder). On the other hand, taking the initial values of the 2-band Cooper 
and forward interactions given by the antiferromagnetic phase in Eq. (B.8), 
one can write down an exact low-energy BCS Hamiltonian of the form [36]: 

H = Ho + Y.J dxVikAlAk, (64) 

i,k 

where V^fe < is consistent with the symmetry is dx2-y2 and the summation 
is restricted to bands in the nodal region, i.e., the umklapp scattering and 
antiferromagnetic processes for those bands have been cutoff at low energy by 
the finite chemical potential. In the large A^-limit phase coherence between 
the bands i and k is due to a Kohn-Luttinger type attraction [37] mediated by 
the short-range antiferromagnetic fluctuations [38,39,40] (In the small A^ limit, 
the phase coherence must be obtained by introducing by hand a Josephson 
type coupling between band pairs). Since the antiferromagnetic processes are 
processes involving 4 bands, this suggests that the antinodal regions should 
contribute to the emergence of superconductivity in the nodal regions. In Sec. 
4.4, we invoke a microscopic Andreev-like scenario for the emergence of d-wave 
superconductivity in the nodal regions mediated by the antinodal regions. 
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Finally, close to "optimal doping", i.e., when there is phase coherence between 
all the band pairs, bosonization predicts [36] 

(Al(x)Afc(O)) oc x-V^, (65) 

as a result of the interband phase coherence. 

4-3 Very large doping: Reminiscence of the Landau- Fermi liquid 



For large enough dopings, all the 4-band and umklapp interactions can be 
neglected. Without umklapp and antiferromagnetic interactions, the weak- 
coupling fixed point of the large N limit is reminiscent of the Fermi liquid. 

More precisely, Lin et al. [88] have discussed the RG flow of the A'^-leg Hubbard 
ladder when including Cooper and forward interactions only. A band pair can 
still develop a spin gap and phase coherence; however, there is no phase coher- 
ence between the band pairs and as a result the gap vanishes exponentially as 
a function of N. More precisely, the forward scattering now only have a small 
phase space in two dimensions [87] and thus can be neglected. In the limit of 
N — > -l-oo, the 2-band Cooper processes give (4c?- = cfj — cij): 

This is reminiscent of a two-dimensional Fermi liquid where one gets: 
dV {91,02) 1 



dl vp 



J dev{ei,e)v{9,92), (67) 



and the angle 9 parametrizes the 2D Fermi surface. This is a renormaliza- 
tion equation for a function, rather than for a constant, i.e., one here has an 
example of "functional renormalization group" . 

This emphasizes that an RG instability can only occur at some points in the 
phase space where one gets attractive or antiferromagnetic interactions. 

It should be noted that for large doping levels, the weak-coupling fixed point 
of the quasi- ID theory is a Fermi liquid because the Luttinger exponents in 
the charge sector converge to unity when the Cooper interactions flow to zero. 



4.4 RG in two dimensions, phenomenology of the pseudogap phase, two gaps 



A major problem in two dimensions is the infinite number of interactions. 
On the other hand, Fermi surface patching [89] is a method to study RG 
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flows of fermionic systems in dimension larger than one. First, we will review 
the 2-patch model which is appropriate when the Fermi surface touches the 
saddle points at (vr, 0) and (0, vr). More precisely, when including a next nearest 
neighbor hopping —t' such that t' > 0, this modifies the band structure as: 

e{k) = — 2t(cos kx + cos ky) + 4t' cos k^ cos ky. (68) 

For if:' = 0, the Fermi surface surface at half-filling corresponds to the umklapp 
surface drawn in Fig. 9. The main effect of t' is to curve the Fermi surface. 
Therefore, when doping the system with holes, one can reach a situation where 
the Fermi surface touches the antinodal (or saddle) points. In the cuprates, one 
can estimate t'/t ^ 1/4 {t' is supposed to vary somewhat from compound to 
compound [90]). For this situation, the Fermi energy intersects the van Hove 
singularity at a hole density ~ 0.2. Moreover, the leading (van Hove type) 
singularity now only arises at the saddle points when adding a next nearest- 
neighbor hopping —t' and t' > 0. On can then simplify the RG analysis by 
including only the antinodal points, resulting in the 2-patch model. An elastic 
umklapp process, where two electrons initially near the saddle point at (vr, 0) 
are transferred to the other saddle point (0, vr), can take place. 

Although it is difficult to build a low-energy theory, one can focus on the 
leading divergent susceptibilities to predict the low-energy phases; one can 
also use numerical techniques as shown by Lauchli et al. [46]. Here, we begin 
to review the results including the effect of the two-particle umklapp scattering 
in the two-dimensional Hubbard model close to half-filling as well as effects of 
the van Hove singularity [41,42,43,45], in the 2-patch model. 

In two dimensions, at the saddle points (see Fig. 10), the susceptibility for 
the (particle-particle) Cooper channel at g = shows a log-square divergence 
y^^-Q ~ ^ (instead of In lu in one dimension) at the energy lu due to the van 

Hove type singularity. For the (particle- hole) Peierls channel a.t Q = {tt,tt), 
there exists a crossover: x^^^ oc In^tu when u > t' and x^^i^) oc Incj when 



uj <^ t \ I . The Peierls channels at g = and the Cooper channel aX q = Q 



also diverge log-linearly but have smaller coefficients and can be neglected. 

Similar to the 2-leg Hubbard ladder, one can identify the forward g2 and (74, the 
backward gi, and the 2-particle umklapp g^, processes. The relevant processes 
are illustrated in Fig. 10. 



The particle-particle Lindhard function at ^ = is unaffected by t': one still gets 
e{k) = e{—k). On the other hand, the particle-hole Lindhard function at g = Q is 
affected at low energy because the t' term results in e{k + Q) / — e(A;). Close to half- 
filling, this will allow the emergence of a new phase dominated by umklapp processes 
(consult Fig. 10) and characterized by the absence of long-range antiferromagnetism, 
similar to the two-leg Hubbard ladder at half-filling. 
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Fig. 10. Scattering processes in the 2-patch model when the Fermi surface hes at 
the saddle points (for t'/t ~ 1/4, this approximately happens for 5 ~ 0.2 where 
< di{y) < 1): forward 54 and 52, backward gi, and umklapp process 53. In the 
half-filled and lightly doped case, Cooper, umklapp, and antiferromagnetic interac- 
tions are strongly coupled. Antiferromagnetic processes are shown in Fig. 9. 

The RG equations take the form (we use the conventions of Refs. [43,45]): 
gi = 2di{y)gi{g2- gi) (69) 

g2 = di{y) (gl + gl) 

gs = -25'35'4 + 2di{y)g-ii2g2 - gi) 
gA = - {gl + gfj ■ 

Here, gi = dgi/dy where y = ln^(a;/£'o) oc 'X^I'_^{lo) and Eq is a high-energy 
cutoff. The function di{y) describes the relative weight between Cooper and 
particle- hole channels: di{y) = dx^/dy. The asymptotic forms are given by 

di{y) 1 when y ~ 1 and di{y) — > ln{t/t') / ^/y when y ocP^ [45]. 

Similar to the 2-leg ladder [26], the forward (inter-patch) coupling strongly 
renormalizes the umklapp as well as the Cooper channel. 

The case di = 1 arises in the limit t' = or in the strong-f/ limit where t' is 
irrelevant in the Peierls channel divergence. This case has been analyzed by 
Schulz [41] and by Dzyaloshinskii [42]. At half-filling, the most divergent term 



At small u), one gets oc Intj whereas x^^^ oc In^w. This ensures di{y) oc l/-y/y. 
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is the antiferromagnetic susceptibility (which has the same exponent as d-wave 
superconductivity but is dominant due to the next leading divergent terms). 
Progress in this direction has been done through the functional RG [44]. It is 
relevant to stress that in two dimensions, it is difficult to keep all the diverging 
channels and build an exact low-energy theory (similar to the quasi-lD case). 
Nevertheless, one can use numerical approaches [46]. It is also interesting to 
observe that the umklapp scattering —>■ +00 strongly diverges, showing the 
occurrence of a Mott state at arbitrarily small U in two (and one) dimensions. 

For finite doping, the flow favors the d-wave superconducting state and the 
antiferromagnetism is progressively weakened [41] . More precisely, at a critical 
doping that depends on the interaction strength, the growth in the antiferro- 
magnetic channels gets cutoff at low energy and the d-wave pairing becomes 
the dominant instability. The "crossover" from antiferromagnetism to d-wave 
superconductivity, for t' = 0, has been analyzed using a N-patch model [91]. 

The extreme limit di = (or t' — > t) has been analyzed by Dzyaloshinskii [42]. 
The RG flow takes the simple form: 



It is relevant to redefine new variables g± = g^^g^ such that the RG flow takes 
the simple form: g± = —g'^. The initial conditions are given hj g^ = g^ = U 
or g^ = 2U and g^ = 0. At low energy, we infer that g-{l) = and (?+(/) 0. 
The extreme limit di = suggests the possibility of a weak-coupling in two 
dimensions close to half-ffiling. On the other hand, one can easily show that 
this weak-coupling fixed point is unstable towards a small but finite value of 
di. More precisely, g- can now access negative regions that allows a strong- 
coupling fixed point governed by g^ —00 and g^ 0. Hence one has 



(73 — i> +00 (and as a result g2 +00) and g4 = —00 at a value yc ~ 

Although one cannot solve for the strong-coupling fixed point using the 1- 
loop RG equations, one can examine the main susceptibilities [45] . It is useful 
to observe that the exponent for the charge compressibility changes its sign 
at di{yc) ~ 0.6. Close to half-filling, this suggests that there exists a critical 
interaction Uc such that for U > Uc the charge compressibility is suppressed 
to zero [45]. Using the results above, one finds Uc/t oc ln~^(t/t'). Close to half- 
filling, this implies a transition from a superconducting phase at U < Uc to a 
phase with a charge gap at U > Uc that can be seen as a precursor of the Mott 
transition. This has also been found in the patch model [44] . It is relevant to 
observe that the fixed point at U > Uc resembles that of the half-filled two-leg 



One gets the precise prefactors via a numerical integration of the RG flows. 



5-3 = -2fi'35'4 



(70) 
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ladder which has a spin and charge gap [46]. A similar phase also occurs in the 
3- leg ladder away from half- filling when including the t' hopping term [34] . 

Note that in two dimensions, starting from the weak-coupling regime, this 
new phase for U > Uc only occurs as a result of the finite value of t'. In two 
dimensions, t' is an important parameter that allows to curve the Fermi surface 
and therefore produces umklapp scattering only in certain directions of the 
Fermi surface. In contrast, in the quasi-lD band model, the band structure 
reveals a hierarchy of velocities, allowing for a truncated Fermi surface close 
to half-fiUing even when t' — 0. At a, general level, it is interesting to observe 
that new spin phases may occur as a result of "frustration" [92] {i.e., as a 
result of the next- nearest neighbor hopping in the weak U limit) . 

It is not so obvious that one can describe the resulting state by a simple 
mean field or Hartree-Fock factorization. In the case of the SO(N) Gross- 
Neveu model arising at the fixed point of the weak-coupling Hubbard ladder, 
a mean-field description is only exact in the large even N limit; for finite N, one 
has to take into account that the gap fiuctuates leading to new bound-state 
excitations, for example. Moreover, it is very difficult within this RG approach 
to clarify how the charge and spin gaps spread out laterally along the umklapp 
surface. Based on the experimental results of high-Tc cuprates [15,16] and our 
results on the N-leg ladder (see Sec. 2 and 3) [30], it sounds legitimate to 
argue that by increasing the electron density (remember that for t'/t = 1/4 
the Fermi surface intersects with the saddle points for 6 = 0.2, close to optimal 
doping) the spin and charge gaps will propagate laterally along the umklapp 
surface that will let a Fermi surface consisting of 4 Fermi arcs centered at the 
nodal points (±7r/2, ±7r/2). This is consistent with high temperature series 
by Puttika et al. [93]. This is also consistent with the Luttinger theorem, that 
tells us that the area enclosed by the arcs inside the half-filled surface gives 
the hole density (see Sec. 4.6). A truncation of the Fermi surface has also 
been found within the SU(2) slave-boson approach of the t — J model (consult 
Appendix C and Refs. [11,13,14]). 

Within this scenario, the pseudogap phase is characterized by (heavy) pre- 
formed pairs and superconductivity might emerge from the proximity effect 
between the Landau quasiparticles of the Fermi arcs and the preformed pairs. 
The pseudogap is not the superconducting gap. Due to the vanishing of the 
charge compressibility in the antinodal regions, the superfluid weight will be 
that of the Fermi arcs alone. The RG approach in two dimensions suggests the 
existence of two distinct energy scales (or gaps) close to the Mott state: T* 
is the temperature at which umklapp scatterings produce a charge and spin 
gap in the antinodal regions, resulting in an (insulating) RVB state, and at 
a lower temperature Tc the Fermi arcs will mediate d-wave superconductiv- 
ity. The microscopic mechanism for the occurrence of superconductivity may 
be an Andreev process (proximity effect), involving the antinodal (hot) re- 
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gions [82,94]: two electrons from the nodal region might virtually 'hop' into 
the antinodal region and may be converted into a Cooper pair as a result of 
short-range RVB correlations in the antinodal region. In the pscudogap phase, 
the specific heat is expected to vanish hnearly as 5 — > (only the Fermi arcs 
contribute) , in agreement with experiment [95] . One important consequence of 
this theory is that the energy gap in the nodal regions does not scale with the 
pseudogap at (vr, 0). In fact, this is in agreement with angle resolved photoe- 
mission spectroscopy (ARPES) measurements [96,97]. Other support for two 
gaps comes from Andreev refiection studies [98] and Raman scattering [99]. 

At this point, it should be noted that the slave-boson approach applied to the 
t — J model gives a distinct scenario for the pseudogap phase (see Appendix 
C). The large gap at (tt, 0) is a d-wave superconducting gap and the pseudogap 
corresponds to a superconductor destroyed by phase fiuctuations. At Tc, the 
holes (holons) condense producing a d-wave superconductor; in this theory, 
that relies on spin-charge separation, there is a single d-wave gap. 

A theory of Tc for the underdoped cuprates has been built by analogy to 
the Kosterlitz-Thouless transition in two dimensions. More precisely, in 1995, 
Emery and Kivelson proposed a model based on phase fiuctuations [100]. They 
pointed out that Tc of the underdoped cuprates is determined by the energy 
scale for phase fiuctuations. In their picture, the entire region below T* is 
characterized by phase fluctuations and at Tc, the phase stiffness Kg jumps 
between zero and a finite value Ks{T~) given by the universal relation: 



Here, m represents the (effective) mass of an electron or 2m embodies the mass 
of a Cooper pair. The precise value of Tc thus depends on the superfluid density 
Ps- Empirically, when the concentration of holes in the Cu02 planes is not too 
large one expects ps{T = 0) oc 5. In fact, this relation is automatically satisfied 
if the Fermi arcs only become superconducting below Tc (see discussion above) . 
This also happens in the renormalized mean-field theory of the t — J model 
presented below. This suggests that Tc oc 5 at low doping levels, as confirmed 
experimentally [101]. On the other hand, as mentioned by Lee and Wen in 
1997 [102], the thermal excitation of quasiparticles near the nodes produces 
a linear decrease of ps with temperature, wich was the earliest experimental 
evidence for d-wave symmetry [103]. The number of quasiparticles in the d- 
wave state varies as since the density of states is linear with energy. On 
the other hand, the amount of decrease of Ps per quasiparticle is inversely 
proportional to its energy, canceling one factor T: thus, Ps decreases linearly 
with T. Following Lee and Wen [102] and MiUis et al. [104], one predicts: 



ksTc = -K,{T-) 



ps 



(71) 



8 m 



p,(T) _p,(T = 0) 
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(72) 
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and the parameter a is given by (see Appendix C.4): 



2 In 2 o'^F 



(73) 



a — 



a — . 

TT V± 



Here, v± represents the velocity of the nodal quasiparticles in the direction 
of the maximum gap. The ratio vf/v± can be measured through the thermal 
conductivity [14,105]. The only assumption made is that the quasiparticles 
with a wave- vector k carry an electrical current j{k) = —eavp, and a is a 
phenomenological Landau parameter [104], which was left out in the original 
paper by Lcc and Wen. Moreover, Vp/v^ is known to go to a constant for 
small doping 5 and ~ 0.5 [14]. This is a strong result because it states 
that despite the proximity of the Mott state the nodal quasiparticles carry a 
current which is similar to that of the tight-binding Fermi-liquid theory. 

Therefore, at T = decreases proportional to doping, yet its rate of de- 
crease with temperature does not vanish with doping, but in fact remains 
relatively constant [106]. The insensitivity of the linear T slope in ps to dop- 
ing was demonstrated experimentally [107]. The decrease of Ps to zero, at a 
given doping level, is also sometimes considered to determine [14,102]. Now, 
coming back to the Kosterlitz-Thouless scenario, one must recognize that the 
Ps which controls the transition is not Ps(T = 0) but Ps{T) which is greatly re- 
duced by quasiparticle excitations. By combining the two effects, the decrease 
of Ps{T) becomes faster than linear and eventually becomes infinitely steep. 
But this happens only very close to T^, therefore, the quasiparticle mechanism 
may give a good (better) estimate of [102]. The Lee- Wen mechanism of 
Tc is relevant for the underdoped side of the phase diagram where it offers 
a natural explanation for Tc ^ Ps{T = 0) 6 and holds all the way up to 
optimality. 



4-5 Gutzwiller projected d-wave superconducting state and mean-field theory 

Now, we analyze the strong- interaction t — J model. First, we closely follow the 
paper by Baskaran, Zou, and Anderson [108], and describe a (renormalized) 
mean-field approach that allows to predict the existence of a superconducting 
dome in a quite natural way and give an idea of the pseudogap [109,110]. 

Assuming strong interactions, the t — J Hamiltonian reads: 



H^-t did. 



f^dlgdis + J ^ Si ■ Sj 



(74) 



{r,j) 



where J — At'^/U and S^' — (i^-fd^. 



, and Hi, 



'IS 



dlgdis is the number occupa- 
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tion per spin. For simplicity, we set the hopping term t' = C. 
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To study the ground state and the excited states of the t — J Hamihonian one 
can use a BCS trial wavefunction P\BCS) [2] where 

\BCS) = + v^,dl/_^^^)\Vac); (75) 

k 

P = Yiii^ ~ ^ii^ii) is the Gutzwiller projection operator [67] that eliminates 
all configurations with double occupancy, as appropriate for U +oo, \Vac) 
is the vacuum state, and vj: are the variational parameters satisfying the 
normalization condition + = 1. At half-filling, the energy of the 
projected d-wave BCS state is very close to the ground state energy; the d-wave 
BCS state is in general a better trial wavefunction than the s-wave BCS wave 
function [112]. The projected BCS wavefunction is a natural generalization of 
the BCS state to strongly correlated systems. This projected wavefunction is 
appropriate to describe three phases: a resonating valence bond insulator, a 
Fermi liquid metal (a state with Uj:v'^ = 0), and a (d-wave) superconductor. 

In general, the Gutzwiller projection is cumbersome to implement (but prop- 
erties of the trial wave function can be evaluated using Monte Carlo sampling) 
and a simple approximate scheme was proposed, called the Gutzwiller approx- 
imation [12]. The essential step is to construct an effective Hamiltonian: 

Heff = -tgt dl^dis - fJ^dl^is + Jgs Y ' C^^) 

such that the pojection operator is eliminated in favor of the reduction factor 
gt = 26/{l + 6) in the kinetic term [113]. In the limit of small hole doping, gt 
26; the electron number per site is n = 1 — 6. Additionally, one gets that the 
projection operator enhances spin-spin correlations: gg = 4/(1 -|- 5^) [12]. The 
renormalization factors gt and gg are explicitly determined by the ratios of the 
probabilities of the corresponding physical processes {e.g., hopping processes) 
in the projected state P\BCS) and non-projected state \BCS). 

By analogy to the slave-boson theory [14], the mean- field Hamiltonian obeys: 
Hmf = -Yi^9t+Xij)dlsdjs-J2f^(^isdis + Yl ^ij (4T4i ~ dl^d]^^+h.c., (77) 

where Xij = X = {■^9sJ / ^){d\gdjs) for nearest neighbors and zero otherwise, 
and for the d-wave BCS case, the variational gap obeys 

= A(cos kx — cos ky). (78) 

A more realistic model for the cuprates gives a similar result, as it must be, since 
details of the tight-binding Hamiltonian should not matter at large U [111]. 
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0.2 
Doping 5 

Fig. 11. The d-wave gap A and superconducting order parameter Asc (~ Tc), from 
the renormalized mean-field theory for J/t = 0.2 [12], as a function of hole doping 6 
(close to half-filling the d-wave gap A converges to J; the gap is in unit of 3gs J/ 4) . 



In this section, to simplify the discussion, we assume that A is real. It should 
be noted that the connection between antiferromagnetic fluctuations and su- 
perconductivity becomes flagrant through the mean-field approach since the 
Heisenberg exchange J gives rise to a pairing contribution. 

Here, the BCS theory is governed by the coupled equations [11,12]: 

^k = l9sJT.^>:_U'^, (79) 

t/ k' 



and 



Xk = -\9sJY.lk-k'^^ (80) 



k' 

where = gtej^ - - and (here) = -t7g where 7^ = 2 (cos + cos ky) is 
the Fourier transform of the exchange interaction (simply the nearest neighbor 
result). = J^i + A| has the same form as in the BCS theory. 



These gap equations can be solved [12]. An explicit derivation of the solutions 
is given in Ref. [12]. Interestingly, the d-wave (spin) gap A falls almost linearly 
with the hole doping from a number of order J, and vanishes around 6 = 0.3 
for J/t = 0.2. This quantity may represent the pseudogap, which is known 
to vary experimentally in this way. A more realistic band structure for the 
cuprates (including the next-nearest neighbor hopping t') gives a similar result 
[111]. The renormalized mean- field theory supports the RVB approach, i.e., 
the appearance of a large energy scale representing the spin gap. The spin gap 
was first observed by Nuclear Magnetic Resonance (NMR) at the end of the 



39 



eighties [9]. Its significance was only slowly recognized by the mid-nineties. 

A represents the d-wave gap or the "pseudogap" . On the other hand, one can 
define a superconducting order parameter Asc which obeys [11,111] 

Ala - {BCS\P{4^dld,+i^d,+i^)P\BCS), (81) 

for large distance I; the real superconducting order parameter is the off diago- 
nal long range order eigenvalue of the density matrix. Using the projected BCS 
state, this averaged value is renormalized by a factor gf, implying that the su- 
perconducting order parameter vanishes linearly at 5 = 0. In fact, Asc 6 
at low doping is a general property of projected superconducting wavefunc- 
tions. The local fixed number constraint imposed by P at 5 = leads to large 
quantum phase fluctuations that fatally destroy the superconducting order at 
very low doping levels. The state at 5 = is an insulator with a vanishing 
Drude weight. Remarkably, A^^; bears a striking resemblance to the variation 
of Tc with doping, and was suggested to be a measure a within this ap- 
proach. Eq. (81) suggests that Tc (only) characterizes the emergence of the 
off-diagonal long-range order. However, as pointed out in Ref. [12], below 
the superconducting order parameter Asc could also be defined as: 

Asc = {BCS\Pid,^d,^)P\BCS), (82) 

which is supposed to renormalize with gt (when taking the averaged value using 
the projected BCS state). This definition would suggest a two-gap scenario 
where the superconducting gap is distinct from the pseudogap (in relation to 
Sec. 4.4). The physical amphtude of the superconducting gap A^^ is shown in 
Fig. 11 [12] and grows linearly with S at low dopings. This is relevant to note 
that those predictions regarding the d-wave nature of the superconducting gap, 
made in 1988, came before the experimental confirmations in 1993-94 [5,6,103]. 
The d-wave pairing symmetry was also predicted by the "spin fiuctuation 
theory" (see introduction), based on a more orthodox structure [40,114]. On 
the other hand, the renormalized mean-field theory is different from the spin- 
fiuctuation theory since the latter is based on a Fermi liquid theory and is 
unable to deal with the unusual properties close to a Mott insulator. 

Variational Monte-Carlo methods for projected d-wave states support these 
results [111]. For more details and references, sec the recent review by Ogata 
and Fukuyama [13]. They also provide a quantitative description on the nodal 
quasiparticles ("nodons") in the BCS state. The Gutzwillcr projected d-wave 
superconducting ground state displays sharp nodal quasiparticle excitations 
and is in good agreement with experiment (below Tc) [115]. These low-lying 
excitations are expected to dominate the low temperature thermodynamics, 
transport, and response functions in the superconducting state. 

Within this scheme, the nodal quasiparticles have a coherent spectral weight 



40 



Z that goes to zero as gt but whose Fermi velocity is very weakly doping 



dependentLfJ and remains non-zero as5— s>0 [11,111]. also implies that 

the real part of the self-energy 3?eE obeys |(93?eS/(9c(j| ~ 1/5 and \d'^eLjdk\ ~ 
1/(5 [vp 7^ 0). The effects of the renormalization gt on the superfiuid density ps 
and on the Drude weight [116] is a consequence of the RVB based theories. In 
fact, the renormalization gt applies to any term in the Hamiltonian which is a 
one-electron energy. In particular, the Green's function of the quasiparticles in 
the superconducting state contains a sharp coherence peak at the quasiparticle 
energy on top of a very broad incoherent spectrum and ARPES experiments 
have confirmed that the amplitude of that peak is proportional to 25 [117,118]. 

To summarize briefly, the Gutzwiller projected d-wave wavefunction confirm 
the occurrence of d-wave superconductivity in the strong-coupling t — J model 
for < 5 < 5c where 5c ~ 0.3 for J/t = 0.2; results in the superconducting 
state are in very good agreement with experimental results (for a more detailed 
comparison between theory and experiment, consult Ref. [11]). 



4-6 Pseudogap phase, generalized Luttinger sum rule and electron propagator 



However, the situation from the renormalized mean-field theory is less clear 
in the pseudogap phase. Phase fluctuations have to be considered thoroughly. 

Moreover, there is another difficulty: the proliferation of nearby alternative 
states of different symmetry. One important issue is the evolution to the anti- 
ferromagnetic state at very low doping. On general principles mesoscopically 
inhomogeneous states (stripes) are likely to be stable at low doping on some 
scale. Charge and spin stripe order have been seen experimentally in a few spe- 
cial cuprate compounds, specifically La2_5Ba5Cu04 and Lai6-5Ndo.4Sr5Cu04. 
The highest stripe ordering temperatures occur at 5 = 1/8, where Tc is 
strongly depressed. It has been shown experimentally that the dominant im- 
pact of the (spin) stripe ordering is to electronically decouple the Cu02 planes 
[119]. The data on anisotropic transport and magnetization properties of 
La2_5Ba5Cu04 at 5 = 1/8 report a crossover to a state of 2D fluctuating su- 
perconductivity, which eventually reaches a 2D superconducting state below 
a Kosterlitz-Thouless transition. Interestingly, it appears that the stripe order 
in La2-5Ba5Cu04 frustrates the three-dimensional phase order, but is fully 
compatible with 2D superconductivity. The layer decoupling mechanism has 
been theoretically addressed in Ref. [120]. Additionally, the nature of stripes 
in a generalized t — J model has been recently investigated in Ref. [121]. 

A second class of competing states has its origin in the SU(2) gauge symmetry 
emerging at half-filling, as explained by Affleck et al. [122]. In the undoped 

This stems from the fact that x does not depend much on doping [14]. 
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state, with exactly one electron per site, the presence of an up spin is equivalent 
to the absence of a down spin and vice-versa, allowing independent SU(2) 
rotations at each site. This implies that an undoped RVB state or any Mott 
state can be represented by an enormous number of wave functions before 
projection; other candidates are the staggered flux state [4] and the d-density 
wave [123]. Wen and Lee [51] have proposed that SU(2) rotations, connecting 
fluctuations of staggered flux states and d-wave superconductivity, may play a 
crucial role in explaining the pseudogap phase (see Appendix C). Very recently, 
Kaul et al. have discussed the possibility of an algebraic charge liquid (in this 
context, the algebraic correlations reside in the charge sector) which seems to 
arise between the antiferromagnet and the superconductor [124]. 

Another approach to capture the pseudogap region is to start with the iV-leg 
ladder close to half-filling [30] or similarly a doped spin liquid comprising an 
array of two-leg ladders close to half-fiUing, in the weak-couphng regime [50] . 

First, using the exact SO (8) low-energy theory of the weak-coupling regime 
(see Sec. 2.3), one can compute the electron Green's function around each 
of the four Fermi points of the two-band model. The single-particle Green's 
function at half- filling (/x = 0) is governed by the coherent part [80] : 

<^a(/C, UJ) - ^2 _ £;2(^) _ + ^^'^^ 

where Ea{k) is the bare dispersion in the corresponding bonding or antibond- 
ing band close to the Fermi wave-vector kpa (close to half-filling, the two bands 
are characterized by the same Fermi velocity), resulting in 

E,{k)=Va{k-kFa)- (84) 

The quasiparticle weight Za ^ I and A is the D-Mott single-particle gap [80]. 
While the coherent part has a form similar to the diagonal Green's function 
in BCS theory, there is no off-diagonal component of the Green's function in 
this case since the system is in a Mott insulating state at half-fllling. 

It should also be noted that Ga{k,0) changes sign from positive to negative 
through a zero as k crosses kpa'- this occurs when Ea{k) = 0. This shows that 
the area characterized by 3?eGa(/c, 0) > is unchanged by the interactions. 
This satisfies the Luttinger sum rule (LSR) [125] for the 2-band system: 

2V / dfc = 4VA;Fa = 47rn. (85) 

The Luttinger sum rule, of course, works for a single-chain Luttinger liquid 
[126]. The key point about the LSR as pointed out in the famous textbook 
by Abrikosov, Gorkov, and Dzyaloshinski [127] and more recently by Tsvelik 
[128], is that the sign change from positive to negative values of ^eG{k, a; = 0) 
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Fig. 12. A sketch of the electronic structure in the pseudogap phase by analogy 
with the A^-leg ladder weak-coupling investigation. In the region S the Fermi sur- 
face shown as "heavy lines" is gapped due to umklapp scattering processes. From 
the quasi-lD analysis with t± ~ t, one can rigorously prove that the D-Mott gap 
propagates from the saddle points along the 2D umklapp surface. The generalized 
LSR for the A^-leg ladder shows that the area enclosed by the Fermi arcs (region A) 
and the umklapp surface (dashed lines) should be proportional to the hole density. 

is not restricted to an infinity in ^eG{k, 0) as in the case of a Fermi surface 
of a Landau-Fermi liquid. It can occur through a zero in ^eG{k, a; = 0) as in 
the BCS theory of superconductivity or as in the 2- leg Hubbard ladder [80]. 

In fact, we can naturally extend this analysis to the case of a A^-leg ladder 
close to half-filling which unambiguously demonstrates a possible truncation 
of the Fermi surface; consult Sec. 3 and Sec. 4.2. When doping, the holes enter 
the band pairs successively from the diagonals of the Fermi surface. For a 
small hole doping 6, the insulating antinodal directions yield exactly the same 
fixed point as at half-filling. This implies that the electron Green's function 
in the antinodal regions is still governed by a self-energy of the form [80]: 

^a{k,u)= (86) 

and the gap Aqq is related to the energy scale defined in Sec. 3, corre- 
sponding to the energy scale at which a band pair (a, a) becomes successively 
frozen out and forms a D-Mott state. It is relevant to mention that this form 
of self-energy also emerges in a doped spin liquid consisting of 2-leg Hubbard 
ladders coupled by a long-range interladder hopping [50]. 

Now, it is interesting to discuss the LSR through this large A^-limit. At half- 
filling, the LSR for the A^-leg ladder system is automatically satisfied: 



2 y / dk = AykFa = 2TmN. (87) 

a=l^..,7V^3?eG„(/c,0)>0 ^ 

The Luttinger contour of zeroes in ^eGa{k,0) corresponding to Ea{k) = 
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coincides with the 2D umklapp surface or with the antiferromagnetic Brillouin 
zone, assuming that t±_ = t. Close to half-fiUing, the bands located around the 
antinodal directions remain insulating and as a result have a self-energy of 
the form of Eq. (86). The remaining Fermi surface then can be viewed as four 
arcs (by analogy to the 2D situation): applying the general formulation of 
the Luttinger theorem requires that the area enclosed by the arcs inside the 
umklapp surface gives the hole density, as illustrated in Fig. 12. Even though 
the Luttinger theorem is applicable in various cases and is general, it can also 
break down in specific two-orbital Mott insulator systems [129]. 

By analogy to the ladder systems [30] and to the array of coupled ladders 
[50], taking the parameters from the rcnormalized mean-field theory in Sec. 
4.5, Yang et al. [48] have suggested a similar single-particle Green's function 
in the pseudogap phase of the underdoped cuprates for the large U limit: 

G^^^(^, u) = ^ ^ + G,nc, (88) 

uj-m-\^u\v{^+uk)) 

where [48] 



^o(^) = -2t{5) (cos + cos ky) (89) 

^ (k) = ^o(^) - 4t'(5) cos kx cos ky - 2t"{6) (cos 2kx + cos 2ky) - /j, 
Ajj = A{6) (cos kx — cos ky) . 

(We have noted u = uj+i6). The factor ^(k) —^o{k) has been adapted from the 
t±{k±) term emerging in the electron Green's function for an array of coupled 
ladders [50]. From the renormalized mean- field theory of Sec. 4.5 [48]: 



m^gtt + ^gsJx (90) 

t'i5)^-gtt' 
t"{S)^gtt". 

where t' > within our conventions; see Eq. (68). Here x is a dimensionless 
parameter (we have slightly changed the notations compared to Sec. 4.5). We 
also introduce a 3^*^ nearest neighbor hopping t". The chemical potential /i 
is determined from the LSR on the electron density. It should be noted that 
the only difference with the weak-coupling Green's function from the quasi- ID 
analysis resides in the factor gt of the renormalized mean-field theory. 

In the limit of zero doping, 5—^0, then Qt and the quasiparticle dispersion 
reduces to the spin dispersion and the quasiparticles has a vanishing weight in 
the single particle Green's function. At small but finite 6, similar to the D-Mott 
phase, the zero frequency Green's function G^^^{k,Q) that enters the LSR 
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(a):x=0.05 



(b): x=0.10 



ic):x=0.14 




Fig. 13. From Yang et al. [48]: Contours on which G{k,0) changes sign at various 
hole concentrations 5 are shown in (a)-(e). A summary is shown in (f). The shaded 
area corresponding to ^eG{k, 0) > satisfying the LSR. In the normal pseudogap 
phase, the line connecting the saddle points (tt, 0) and (0, tt) is the Luttinger surface 
of zeroes and the pockets in the thick line represent the infinities of ?R.eG{k, 0). The 
chosen parameters are: x = 0.338, J/t = 1/3, t' /t = 0.3, and t" /t = 0.2. 

has lines of zeroes when ^o(^) = = — 2t(5)(cos /c^ + cosA;^). The Luttinger 
contour of zeroes in G^^^{k, 0) then consists of straight lines connecting the 
saddle points (ivr, 0) and (0, ivr) and perfectly coincides with the 2D umklapp 
surface (see Fig. 12). This Luttinger contour is also in accordance with the 
functional RG calculations [44] (supplied by a numerical treatment at the fixed 
point [46]) on the weak-coupling 2D t — t' — U Hubbard model of Sec. 4.4. 

A second feature that emerges from G^^^{k,Lj) is the appearance of hole 
pockets (rather than arcs) at finite hole doping. The hole pockets define Fermi 
surfaces where G^^^{k,uj) changes sign through infinities and contain a to- 
tal area equal to the hole density. The LSR is satisfied since the area with 
^eG^^^{k, 0) > is bounded by the Luttinger surface which contains exactly 
1 electron per site, minus four hole pockets which have a total area related to 
the hole density. The Fermi surface of the four hole pockets is defined as: 



(91) 



The evolution of the hole pockets as a function of doping is shown in Fig. 
13. The chemical potential /i was adjusted at each 6 to reproduce the correct 
area for the hole pockets. One can see that the hole pocket evolves gradually 
into a more normal surface as S increases. It is also instructive to observe that 
within this ansatz of G^^^ the spectral weight of the quasiparticle pole varies 
strongly around the pocket [48]. In fact, the quasiparticle weight in G^^^ have 
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very small values on the back sides of the hole pockets, which can account for 
the failure to observe these parts from photoemission experiments. However, it 
should be noted that Fermi pockets have been recently inferred from quantum 
oscillation datas [130]. Another explanation based on the algebraic charge 
liquid can be found in Ref. [124]. Many experimental features, e.g., the slow 
variation of the nodal Fermi velocity vp and the Drude weight scaling with 
5 are reproduced by the propagator of Eq. (88). The ansatz for G^^^ agrees 
with ARPES and resolves an appearent disagreement with the LSR [48]. 

The Green's function G^^^ also produces zeroes in the nodal regions along 
the 2D umklapp surface: in fact, the zeroes only disappear at the nodal points 
since the d-wave RVB gap vanishes at those points. The umklapp surface 
along which the energy gap opens up lies above the chemical potential, caus- 
ing particle-hole asymmetry in the quasiparticle spectra in the pseudogap 
state. The spectra are symmetric only along the nodal directions where the 
gap vanishes, and the asymmetry increases away from the nodal directions in 
agreement with ARPES measurements. For a recent discussion, see Ref. [131]. 

Below Tc, following the quasi-lD analysis of Sec. 4.2, the Fermi arcs now 
become superconducting. Therefore, the self-energy in the superconducting 
phase must be modified to convert the infinities of the single particle Green's 
function at the Fermi surface of the four hole pockets onto zeroes, except along 
the Brillouin zone diagonals where the superconducting gap /S.sc{k^) vanishes. 

In the publication of Ref. [48], Yang et al. have proposed the following ansatz 
for the d-wave superconducting state (pushing forward the "2 gap" idea): 

E5c(^,^) = ^rvb{K^) + J^^^(^)l' ^, (92) 

and T^RVBikjUj) is the RVB spin liquid self-energy in Eq. (87): 

^RVsiK^) = \Aj:\y{uj + ^o{k)). (93) 

Gonsult also Ref. [132]. Following Yang et al, the coherent part of the single- 
particle Green's function in the d-wave superconducting state then reads: 

^^-'^^^ ^ — V ^^^^ 

This form of Green's function has been found by solving the coupled equa- 
tions which connect the regular {G^^^^ = — i((i|(i|)) and the anomalous (F^ = 
—i{d^d\^)) Green's functions (following Abrikosov, Gorkov, Dzyaloshinski, page 
295 [127]). In the superconducting state, the off-diagonal component of the 
Green's function becomes non-zero. Including the pseudogap in the normal 
phase and ignoring the effect of the Gutzwiller projector, then one gets [48] : 
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{u - i{k) - Y,nvBM)Glh - iAscF}{k,u;) = 1 (95) 
(<^ + ^(k) + '^RVB {k, -u;))F^Ak,oo) + iAlcGUk,co)^0. 

The superconducting gap function has a d-wave form and is related to the 
anomalous Green's function in a standard way: 

Ascik) = J dk'dwg{k-k')F{k',u), (96) 

where g{k — k') is the d-wave pairing interaction. Within the renormalizcd 
mean- field theory, any single-particle operator must be renormalizcd by gt- 
This justifies Eq. (93) (which includes the self-energy of the pseudogap phase). 

First, since "E^yBik, 0) +oo on the surface where ^o(^) = 0, Geo?* continues 
to have a Luttinger surface of zeroes in the superconducting state on the 
umklapp surface. However, there is an additional set of Luttinger surfaces 
defined by ^{k) + T,RVB{k, 0) = 0. But these are just the Fermi surface of the 
four hole pockets in the normal phase which have now been converted to a 
Luttinger surface of zeroes in the superconducting state. The G^oA satisfies 
the LSR. We also remark that along the Brillouin zone diagonals both and 
Asc{k) vanish, and exactly along these directions there is only a quasiparticle 
pole which crosses the Fermi energy at a Fermi wave- vector given by ^{kp) = 0. 

In a more recent publication, Yang et al. [49] have observed that in the absence 
of particle-hole symmetry in the normal phase, it is more rigorous (justified) 
to use the proper quasiparticle basis in the normal phase such that: 

Gloh = Yl A2 (k) ' ^^^^ 

k 



where = ± + el, e| = ^kMh) + I A,-p, = m - Uk))/2, and 
= 1 + |A^p/(£'^ + ^o(^))^- In the normal phase, the hole pockets 
satisfy the equations E"^ = that is equivalent to e| = or Eq. (91). In 
the superconducting phase, the four hole pockets are naturally converted into 
Luttinger surfaces of zeroes. In fact, this propagator in the superconducting 
phase gives a better agreement with the experimental results and in particular 
with Scanning Tunneling Microscopy (STM) results of Kohsaka et al. [133]. 
Especially, the asymmetry in the tunneling density of states then becomes 
weaker due to the inherent particle-hole symmetry of superconductivity. 

At this point, it is relevant to compare the phenomenological forms for G^^^ 
and Gco^i other proposals. It is very close to the form derived by Ng [134] 
based on spin-charge separation but with an added "not-too -rigorous attrac- 
tion between spinon and holon (see discussion in Appendix C) which leads to 
quasiparticle poles and a self-energy similar to Eq. (97). Earlier, Norman et 
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al. [135] introduced d-wave pairs plus lifetime broadening in the pseudogap 
phase. The effect of classical phase fluctuations on the quasiparticle spectra 
has been addressed by Franz and Millis [136]. A recent comparison between a 
number of theoretical proposals for the arcs has been drawn in Ref. [137]. 



4-7 Overdoped cuprates and Breakdown of Landau- Fermi liquid 

Now, we discuss the crossover to the overdoped regime. Prom the A'"-leg ladder 
analysis of Sec. 4.2 and 4.3, one can conclude that the Fermi hquid is stable 
at low energy when antiferromagnetic processes vanish completely, as a re- 
sult of the circular Fermi surface that forms when overdoping the system. In 
this sense, the onset of superconductivity with doping should coincide in the 
normal phase with the onset of antiferromagnetic (umklapp) scattering. 

Recently, Abdel-Jawad et al. [138] have reported a correlation between charge 
transport and superconductivity in heavily overdoped high-Tc superconduct- 
ing cuprates. For a review, consult Ref. [139]. They found that the onset of 
superconductivity with doping coincides with the appearance in the normal 
phase of strong anisotropic quasiparticle scattering. The application of a mag- 
netic field to suppress the superconductivity revealed that the anisotropic term 
in the inplane (ab) transport scattering rate is linear in temperature violating 
the quadratic dependence characteristic of a Landau-Fermi liquid. 

Through a careful analysis, the authors report a scattering rate of the form: 

r = Fo + aT^ + bT cos2(20), (98) 

where is the azimuthal angle. The isotropic part is standard for metals: the 
sum of impurity scattering (Fq) and electron-electron scattering (oT"^). The 
last term, the anomalous one, is of particular interest. The linear T dependence 
of the in-plane resistivity stems from this term. In Tl-2201 samples [138] where 
doping 6 = 0.25, the resistivity is not perfectly linear, but best described 
by: p — Pq + AT^ + BT [140,141]. The existence of a large Fermi surface 
in the overdoped regime has been seen through the observation of quantum 
oscillations [142] . The key finding of Ref. [138] is that the anomalous scattering 
is deeply anisotropic. It goes to zero as = 7r/4 and is maximum at = vr. 
This angle dependence suggests that the anomalous scattering has the d-wave 
symmetry. It should be noted that the linear T term disappears at S — 0.3 
[140,141]; it coincides roughly with the onset of superconductivity at 5 — 0.27. 

Earlier theoretical investigations of the 2D square lattice via the functional RG 
method found that d-wave pairing in the overdoped region of the phase dia- 
gram was driven by the appearance of a strongly anisotropic scattering vertex 
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Fig. 14. Phase diagram inspired from Taillefer [139]. Three regions can be distin- 
guished from the behavior in the normal phase. For very large dopings, one gets a 
fully cylinrical Fermi surface, no superconductivity, and an isotropic scattering rate 
{5 > 0.27). The overdoped region (0.15 < 6 < 0.27) marks the simultaneous appear- 
ance of anisotropic d-wave superconductivity and anisotropic scattering mediated 
by antiferromagnetic processes at (vr;7r). In the underdoped region, when the van 
Hove singularity touches the saddle points then the Fermi surface will become trun- 
cated in the antinodal directions where the scattering is maximum; the antinodal 
regions will form an insulating RVB region for U >Uc where Uc/t oc ln"^(t/t') (see 
Sec. 4.4). By increasing the electron density, the charge and spin gaps propagate 
along the 2D umklapp surface and satisfy the generalized LSR (see Sec. 3 and 4). 

in the particle-particle and particle-hole channels [44,143,144]. Further inves- 
tigations revealed that the self-energy is also anisotropic [145,146,147,148]. 

Recently, Ossadnik et al. [149] have presented an extensive study of the doping 
and temperature dependence of the quasiparticle scattering rate with the pair- 
ing instability suppressed. Within the RG approach, it is difficult to introduce 
a magnetic field, therefore the d-wave pairing instability has been suppresed 
by an isotropic elastic scattering term added in the free part of the action. 

The self-energy at the Fermi surface is given the by the expression: 

Ea=o = r rfA e{\i{k{)\ - A)5(|e(4)| - A)^(A - (99) 

J ho 

X VlGo{h)Go{k^)Go{h), 

and summation and integration over internal momenta and Matsubara fre- 
quencies is implied. The two-loop diagrams contributing to the self-energy are 
shown in Fig. 6 of Ref. [146]. Here, ki = {ki,LJi), ^{k) = e{k) — where e{k) is 
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given in Eq. (68), and Go is the free propagator. The first propagator contains 
a sharp infrared cutoff xa = ^(l'C(^i)l ~ ^)) fhe second has support at, and 
the third below the cutoff A. The initial condition gives Voo{ki, k-^) = U 
and the frequency dependence of all the vertices has been neglected. To com- 
pute the self-energy, one needs to evaluate the 4-point vertex Va(A;i, /c2, ^3) 
at the energy scale A. As we focus on scattering rates at the Fermi surface 
or 'ismTi\=Q{kF,uj), for a frequency Yimui->Q^ + with a; = 0, one gets an 
imaginary part oc 5{^{k^) — ^(/c2) — reflecting energy conservation [146]. 

It turns out that neglecting the flow of the 4-point vertices, or setting Va. = U , 
is equivalent to a second order perturbative calculation of the scattering rate, 
which gives a behavior away from the van Hove singularities. All deviations 

from the Landau theory scaling form then must be attributed to the renor- 
malization of the 4-point vertices. Now, let us imagine that at energy A, we 
focus on the 4-point vertex Va.(A;i, A;2, /cs), and we fix the outgoing wave- vector 

— * 

close to (tt, 0); Ossadnik et al. have shown that the strongest scattering 
process occurs at (7r,7r). They also predict that these scattering vertices are 
roughly oc A~^/^ for small values of the cutoff A ~ T [149]. Moreover, for very 
small cutoff, the square root divergence is suppressed due to the presence of 
the isotropic scattering rate added in the free part of the action. 

The factor in the self-energy calculation then results in QmE oc T [146]. 
It should be noted that this result is not due to a proximity to the van Hove 
singularity at the saddle point of the band structure: in this case, this lies well 
below the Fermi energy. In contrast, the breakdown of the Landau-Fermi liquid 
theory here emerges from the increase of the 4-point vertex with decreasing the 
energy scale. It should also be noted that this increase in the 4-point vertex is 
not restricted to the Cooper channel; examination of the RG flows shows that 
several channels in the 4-point vertex grow simultaneously, e.g., the particle- 
particle umklapp process. On the other hand, it is quite surprising that this 
increase of the 4-point vertex gives exactly the same profile of $5mS oc T as 
when the van Hove singularity hes at the (antinodal) saddle points [146,150]. 

We argue that the phenomenon associated with the increase of the 4-point 

vertex Va when decreasing the energy scale (thus producing the simultaneous 
reinforcement of several channels) can be seen as a precursory effect of the 
RVB behavior taking place at lower doping levels; see Fig. 14. 

To summarize briefly, Ossadnik et al. [149] have found that the scattering 
rates have a minimum in the nodal direction (0 = 7r/4) and they increase 

towards the antinodal direction. Moreover, the coefficient b in the scattering 
rate is found to increase when decreasing hole doping, whereas the coefficient 
a of the Landau-Fermi liquid does not change much with doping. This is in 
agreement with the experimental results of Abdel-Jawad et al. [138]. 
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5 Cold Atomic Fermions and Hubbard model in optical lattices 

The experimental realization of degenerate Fermi atomic samples (see Ap- 
pendix D.l) have stimulated a new wave of investigations of quantum many- 
body systems [151,152,153,154]. Experimental advances include the observa- 
tion of (molecular) fcrmion-pair condensates [155]. Recently, both repulsive 
and attractive Hubbard models have been successfully implemented in opti- 
cal lattices with fermions (using atoms) [56,57]. Fingerprints of the Mott 
state have been observed for repulsive interactions [57,58]. This suggests that 
d-wave superfluidity on a square lattice might be realized in the future; s-wave 
superfluidity in an optical lattice has already been reached (using ^Li) [156]. 
Superfluidity in optical lattices has been discussed by Hofstetter et al. [157]. 

5.1 Light and Atomic Parameters 

First, we find it instructive to discuss the parameters of those atomic Hubbard 
models, resulting from the trapping of atoms in a periodic potential creating 
by interfering laser beams (for a review, see Ref. [158]). 

Consider an ensemble of fermionic atoms illuminated by several orthogonal 
laser fields (tuned far from atomic resonance). These fields produce a periodic 
potential for atomic motion in two or three dimensions of the form (for more 
details, consult Appendix D.2): 

2(3) 

V{x) = VoY.cos\kxi), (100) 

i=l 

where k is the wave- vector of the hght. The depth of the lattice in each direc- 
tion is determined by the intensity of the corresponding pair of laser beams 
which is easily controlled in an experiment. The potential depth Vo is usu- 
ally expressed in terms of the recoil energy E,,. = h'^k'^/2m. It is interesting 
to observe that for a YAG laser with A = 27r/A; = 1.06/xm and ^Li atoms, 
Er ~ 1.4/xK. When venturing in the community of cold atoms, condensed 
matter physicists who are used to express energy scales in Kelvins (or electron- 
volts) need to remember that in units of frequency: l/iK ~ 20.8kHz. 

The Bloch functions depend on the quasi-momentum and on the band index. 
It is instructive to focus on the corresponding eigenenergies for different lattice 
depths Vq/Ej.. Already for moderate lattice depths of a few recoils, one finds 
that the separation between the lowest lying bands is much larger than their 
extend [159]. One can then focus on the lowest Bloch band and define Wannier 
wavefunctions that are optimally localized at the site Xi. We are interested in 
a situation in which there is roughly one atom per lattice site; such atomic 
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densities generally correspond to free-space Fermi energies on the order of Er 
(note that the kinetic energy and Fermi energy of atoms in the optical lattice 
are not the recoil energy but rather the bandwidth of the Bloch band under 
consideration, which strongly depends on the laser intensity Vq). Then, the 
atoms can tunnel from one site to another. 



To realize the fermionic Hubbard model, one requires that two kinds of atoms 
are present; they can differ by angular momentum or by generalized spin (s = 
1,1). For sufficiently low temperatures, the atoms are confined to the lowest 
Bloch band, and the system can be described by the Hubbard Hamiltonian 
given in Eq. (1). The parameter t (depicted in Fig . 15A) corresponds to the 
tunneling matrix element between adjacent sites 
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t^EriA/V^)eexpi-2e), (101) 
where ^ is precisely given by ^ = (Vo/ErY^^, and the parameter [55] 

U = E^a,fcy^^^ (102) 

characterizes the strength of the onsite interaction. The sign of the scattering 
length tts determines the nature of atomic interactions and can be controlled 
(tuned) experimentally: negative corresponds to attraction between atoms 
{U < 0) and positive corresponds to repulsion {U > 0). 

For example, in the case of ^Li both cases can be realized depending on the 
particular electronic states that are being trapped [152,153]. The hopping 
amplitude t and the on-site interaction U, calculated for the lowest band, are 
plotted as a function of Vo/Er in Fig. 15B (for a cubic lattice). 

The realization of the Hubbard model with ultracold fermionic atoms requires 
3 prerequisites: i) neglect the second band, ii) neglect other interactions besides 
the Hubbard interaction U, and iii) replace the interatomic potential by a 
pseudopotential approximation. Assumption i) is justified as long as the second 
band is not populated (less than two fermions per site, and Vq not too small; 
typically, Vq ^ 2.3Er such that the two bands do not overlap). We also assume 
that the on-site interaction U is smaller than the separation between bands 
Ss- The assumption ii) is justified if other interactions are negligible. The most 
dangerous coupling turns out to be a 'density-assisted' hopping between two 
nearest-neighbor sites [160]. The pseudopotential approximation consists to 



' For a description of the Wannier wave functions and their corresponding tunneling 
amplitudes, consult Ref. [55]. The behavior of the hopping parameter t can be qual- 
itatively understood from a Wentzel-Kramers-Brillouin type argument. The value 
of the Hubbard parameter is obtained by taking the Wannier state localized around 
a bottom well as the Gaussian ground state in the localized oscillator potential. In 
the quantum optics community, the hopping term is often called J. 



52 



B , 




Fig. 15. A) Implementation of the atomic Hubbard model. B) The hopping strength 
t and the Hubbard interaction \J can be tuned via the lattice Vq. C) and D) Different 
lattice geometries can be engineered by resorting to various pairs of lasers and by 
controlling the laser intensities. We have introduced the lattice spacing a = A/2. 



approximate the interaction potential by l^! 



int 



(X — X 



g6{x~x') where g 



Airh ag/m. The Hubbard interaction U is obtained by rewriting the interaction 
Hamiltonian in the basis of Wannier functions. Therefore, the assumption iii) 
is justified when the typical distance between two atoms in a lattice well (which 
is given by the extension of the Wannier wavefunction /) is must be larger than 
the scattering length / 3> a^. It is relevant to observe that for deep lattices 



Vo ^ Er this assumption coincides with the requirement thai] I U <^ Ss and 
for large Vo/Er, boils down to: 



< 



Er 



-1/4 



(103) 



For a deep lattice, the scattering length should not be increased too much. 



It is important to emphasize that optical lattices allow the realization of a 
variety of lattice geometries. The lattice site positions Xi determine the lattice 
geometry. For instance, the arrangement of three (two) pairs of orthogonal 
laser beams leads to a simple cubic (square) lattice. Since the laser setup 

IS One can use I ~ a{Er/Voy/^, U/Er ~ {as/a){Vo/ Er)^/^ , and 6s ~ (ErVo)^^^; we 
have re-introduced the lattice spacing a. We have used the model of an harmonic os- 
cillator with a typical frequency co ~ y/ErVo; from the ground state of the harmonic 
oscillator, one finds / ~ {ErVo)~^/'^ (A ~ 1/k ~ l/y/E^)- 
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is very versatile different lattice geometries can be achieved quite easily. For 
example, three leaser beams propagating at angles 27r/3 with respect to each 
other in the xy plane (and all of them being polarized in the z direction) allow 
to realize a 2D triangular lattice; an additional pair of lasers in the z direction 
can be used to create locahzed lattice sites (see Fig. 15D). The motion of the 
atoms can be restricted to one spatial dimension by large laser intensities. 

The size of an optical lattice is relatively small in comparison to its periodicity 
a with typically a few hundred lattice sites in each dimension only. Also, this 
lattice is out into a slowly varying harmonic potential (trap) which translates 
into a space dependent site offset or a local chemical potential e{x) (see Fig. 
15A). This can lead to the simultaneous coexistence of spatially separated 
regions where alternating superfluid and Mott phases are present [159]. 

Below, we consider that the harmonic potential is smooth enough such that a 
large number of sites are subject to an almost fiat potential. Those sites are 
described by the same filling factor which can be controlled experimentally, 
by either increasing the total number of atoms, by reducing t via an increase 
of the lattice depth or via an increase of the overall harmonic confinement. 

A next important question is the temperature of the many-body atomic fermion 
system. It has been shown that one can determine the temperature of a non- 
interacting 50/50 spin mixture of fermions by measuring the number of doubly 
occupied sites. When U = 0, for zero temperature and deep optical lattices, 
one expects all lattice sites to be occupied by a spin-up and a spin-down 
equally. For finite temperatures, atoms can be thermally excited to higher 
lying lattice sites at the border of the system, thus reducing the number of 
doubly occupied sites. By converting doubly-occupied sites into molecules, it 
has been possible to determine the number of doubly occupied sites versus 
singly occupied sites and obtain an estimate of the temperature for the sys- 
tem [161]. Nowadays, one can reach temperatures ~ O.IT^ where Tp is the 
Fermi energy [ks — 1) oi a non- interacting Fermi gas in an harmonic trap. 

When loading the atoms in the optical lattice, the bandwidth of the Bloch 
band under consideration becomes much smaller than the recoil energy Ej. . 
When the lattice is gradually turned on, adiabatic cooling is expected to take 
place (the entropy of the system cannot fiow and the degree of degeneracy is 
conserved). This leads to a (smaller) temperature of the lattice, say, Tiat- Here, 
the temperature here is not determined by some external bath but just by the 
entropy. The density of states is enhanced considerably as the band shrinks 
(since the one-particle states all fit in a smaller energy window); therefore, 
one expects that Tiat is reduced compared to the initial temperature in the 
trap. Interactions can significantly modify these effects; see below. The Mott 
insulating regime can be reached assuming that U 3> {t^Tiati fj)- Recently, it 
has been confirmed experimentally that for U ~ 24i the double occupancy is 
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strongly reduced to values systematically below 2% (for small atom numbers) 
[57]. This constitutes a fingerprint for the suppression of fluctuations in the 
occupation number. Schneider et al. have recently developed a new method to 
measure the compressibility of the many-body system [58]. The entrance into 
the Mott regime is signaled by a minimum in the global compressibility. 

Concerning the antiferromagnetic ground state, a recent discussion on the 
Neel temperature Tjq can be found in Ref. [162]. It is relevant to observe that 
according to quantum Monte Carlo simulations of the Hubbard model on a 
cube [163], for a given nearest- neighbor hopping matrix element the highest 
Tn ~ t/3 occurs for U ~ 8i, while for a given U the maximal Tat ^ U/20 occurs 
at t ~ 0.15U. Thus, to increase Tn one wants to move to large t (which implies 
a weaker optical lattice, or smaller Vq), and to larger U (which means a larger 
scattering length a^). On the other hand, this moves the system away from 
the regime where it is well- approximated by the single-band Hubbard model; 
see Eq. (103). Mathy and Huse [162] have focused on the precise conditions 
which maximize the Neel temperature in the Mott insulating regime. 

It should be noted that numerical simulations [163] predict T^/J ~ 0.957 
on the cubic lattice (quantum fluctuations reduce T^/J from its mean-field 
value). Hence, T^/t or J/t becomes small (as t/U) in the Mott insulating 
phase. Assuming t ~ Q.IU then temperatures Tiat on the scale of ~ lO^^i^r 
must be reached. Experimentally, it is crucial to reach such energy scales to 
prove the disappearance of the huge spin degeneracy below and at the 
same time to observe d-wave superfluidity since Tc is smaller than J. 

In the experiment [58], initial temperatures in the trap are in the range 
0.15 < T/Tf- < 0.2. At these temperatures, the entropy per particle is larger 
than 2 In 2. When switching on the lattice, a Mott insulator with unit fill- 
ing and In 2 entropy per particle form in the center of the trap even for 
T/Tp = 0.15, for which the average entropy per particle is above 2 In 2. This 
is possible only due to the inhomogeneity of the system, as most entropy is 
carried by the metallic shells at the edges (where the entropy per particle can 
be very large). Therefore, the temperature Ti^t ~ t remains larger than T/v 
(but much smaller than U). The next experimental challenge now consists to 
access temperature scales Tiat of the order of the Anderson's superexchange 
J = At^ /U . Experimentally, one may first cool down the gas (in the absence 
of the optical lattice) down to a temperature where the entropy per particle 
takes a low enough initial value (smaller than In 2); interaction effects in the 
optical lattice then lead to adiabatic cooling mechanisms which could help to 
reach the Neel phase as U/t is increased [164]. A detailed discussion on this 
issue can be found in Refs. [160,165] (but the effect of the trapping potential 
has not been discussed). It seems very useful to think in terms of entropy... . 

Nevertheless, it is relevant to mention that the existence of superexchange 
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antiferromagnetic correlations has been shown in an array of isolated double 
wells [166]. Also, the regime Tiat 3> J is appropriate to observe the so-called 
spin-incoherent Luttinger liquid behavior of ID fermion systems [167,168]. 



5.2 Plaquette models with ultracold fermions and d-wave superfluidity 



In fact, it has been known for some time that the minimal 2D unit that can 
sustain d-wave pairing physics is a single square plaquette [59]. 

Similar to the ladder system, the single plaquette model has been another 
relevant starting point of various approaches to the Hubbard model on a square 
lattice [60,61]. Interestingly, even for purely repulsive interactions, there is a 
range of parameters where two holes tend to bind together on a single plaquette 
rather than to delocalize among different plaquettes. Here, we review the main 
results of the plaquette system along the lines of Ref. [59]. Those plaquette 
systems can be realized with atoms loaded in 2D optical superlattices [62]. 

First, let us consider a plaquette of four sites with four electrons. As soon 
as one switches on the Hubbard interaction t/ 7^ 0, it has been shown that 
the resulting ground state is a singlet with d-wave symmetry. For large U, 
the gap to the excited states is ~ J (the splitting between the lowest lying 
singlet and triplet state: "the spin gap"), whereas at small U the (spin) gap 
varies as U'^ [169]. In fact, all eigenstates of a 4-site Hubbard model can be 
calculated analytically by taking advantage of the symmetries of the model 
and the results have been published by R. Schumann [170]. Assuming that the 
inter-plaquette hopping t' is small enough (see Fig. 16), then the gap can be 
treated as big. When the plaquette contains two electrons, the ground state has 
also been shown to be a singlet as a result of antiferromagnetic correlations. 

Now, let us discuss the emergence of d-wave pairing. For non-zero U, this can 
be seen by writing the wave function with two electrons (on 4 sites) as [59]: 

1^2) = K (4^4^ + dl^d\^ + ...) \Vac), (104) 



where K is a normalization facton I. The (1, 2) amplitude has the same sign 



as the 90° rotated (1,4) amplitude. One can thus write: 

1^2) = Al\Vac), (105) 
where Aj is an operator that creates an s-wave pair. On the other hand, when 



The signs are in agreement with the wavefunction at [/ = 0: [^'2) oc (djj^ + d^^ + 

+ c^4|)(c^J| + (^21 + ^3t + d\^)\Vac). If a nonzero U is added to the Hamiltonian, 
all ampUtudes remain positive, although they no longer have the same amphtude. 
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Fig. 16. A plaquette is the minimal system that exhibits d-wave symmetry. When 
loaded with 4 fermions the ground state is d-wave symmetric whereas with 2 elec- 
trons the ground state exhibits s-wave symmetry. Consequently those two states 
are connected through the d-wave pair creation operator. An atomic checkerboard 
lattice may be realized by coupling the plaquette through an hopping amplitude t' . 

one starts from the Mott state the ground state is not the vacuum but rather 
the Neel state with exactly one particle per site Y^a)- The interesting point 
is that this same two-particle ground state {^2) can also be created by a 
d3;2_j^2-wave operator removing particles from the Mott insulating state. For 
the model system with four sites, (\1/2| A(i|\l/4) is large whereas (\[^2|As|^4) = 0. 

More precisely, for a repulsive U the largest real space amplitudes in the 4- 
particle ground state wavefunction are for the "Neel" configurations: 

|^]v) = 4i4i4T<l^«c), (106) 

and the spin reversed state 

1^2^) = dl^d\^d\^d\^\Vac). (107) 

The emergence of d-wave pairing in the plaquette can be seen as follows. We 
examine the relative phase for an electron pair on sites (1,2) and the 90" 
rotated pair on sites (1,4). Annihilating the appropriate electrons, one gets 

^31^41 1 ^]v) = {d^^di^) {d\/^^dl^d\^) \Vac) = -dl^d\^\Vac), (108) 

and: 

4|rf2i|^]v) = {d^]d2i) {d\^d\^dl^d\^) \ Vac) = +dl^d\.^\Vac). (109) 
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Fig. 17. In principle, different lattice geometries may be realized through superlat- 
tices [171]. (A) Ladder system. (B) Checkerboard lattice around 3/8 filling. 

Thus, to have a non-zero overlap against the state \'^2), one must use 

Ad = {d3^d2i- d3id4i + ...), (110) 

on the four-electron Mott insulating state; because of the minus sign, this is 
a dx2^y2-wave operator. To check that is a singlet one can operate on the 
linear combination |\l/]v) + I^at)) the 4-particle ground state, |\I/jv) I^a^) 
enter with a relative plus sign. The linear combination which gives the same 
relative signs as in |\l/2) is obtained using (the d-wave operator) [59] 

[(4TC?2i - 4irf2T) - {ds^d^ - ds^d^^) + ...] (|*]v) + 1^^))- (111) 

To summarize, the hole-pair creation operator connecting the |\[^4) and \'$2) 
state must have a d-wave symmetry. 

When the hopping between plaquettes is zero, t' = 0, thus the ground state of 
the model is known exactly for all U/t and (also all band- fillings or dopings). 
In fact, the ground state of the single plaquette problem can be obtained either 
by brute force [170] and Bethe Ansatz [169,172] (four site plaquette is a ID 
ring with four sites). At exactly half- filling, the ground state is an insulator 
for U/t > (similar to the 2D Hubbard model and the ladder systems). 

Additionally, the insulator has a d-wave symmetry and is called the D-Mott 
plaquette state [173]. It should be noted that this D-Mott plaquette state is 
slightly different from the D-Mott state for the two-leg Hubbard ladder system: 
here, the plaquette state is odd under 7i/2 rotations about the plaquette center. 

Now, let us look at the hole-doped case, with 2 holes in two isolated plaquettes. 
Similar to the ladder, the holes can bind together within the same plaquette 
or separate as single holes in each plaquette, depending on the binding energy: 

A, = 2Eg{3)-Eg{A)-E,{2), (112) 
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being positive or negative, respectively. Here Eg{Na) is the single plaquette 
ground state energy when loaded with Na atoms. A positive Af, signifies an 
efi:ective attraction between doped holes in the sense that for two doped holes, 
it is energetically preferable to place both on one cluster than to place one 
on each of two clusters. For the isolated square, it has been shown that A;, 
is positive for < ?7 < ?7c ~ 4.58t [170]. Then, two adjacent plaquettes can 
be coupled through a weak tunneling t' to form a super-plaqucttc. As long as 
< t' < Af), the states with 4 particles in one plaquette and 2 particles in the 
other are lower in energy and the states with one hole in each plaquette are 
energetically suppressed. Thus, starting with the half-filled case and doping 
with holes, may result in d-wave superfiuidity since (^'2|Ad|^'4) is large [62]. 

In fact, a quite ideal situation to observe d-wavc superfiuidity and avoid occu- 
pation states of a plaquette with 3 particles is to tune U 2t ( A^ is maximum) 
and to consider situations close to the 3/8 fiUing, where the number of 4 par- 
ticle states and the number of 2 particle states are (almost) equal. 

More precisely, by treating the |\1'2) and |\l/4) states as the pseudo-spin compo- 
nents I t) and I I) of an effective spin system, in the limit where A^ 3> t'^ where 
t'g is the hopping amplitude for each specie, then one can build an effective 
XXZ Hamiltonian (along the lines of the ladder case of Sec. 2.1): 

Heff = E {-M'^n'^M' + 44)) + ^^44" (113) 
where the summation involves two neighboring plaquettes and [62]: 

J^-/-^, J' = itf + tf)^^; (114) 

In fact, the parameters g{U /t) and gz{U/t) have been properly derived in Ref. 
[62]; see Appendix D.3. At 3/8 filling, the behavior of the function gz{U/t) is 
quite complicated. On the other hand, for t'-^ ~ t'^ = t', the system will be in the 
d-wave superfluid phase as long as g ^ gz' this requires <^ U < 2. It. Away 
from the 3/8 filling, the propagation of the hole pairs (plaquette states) in the 
system will be always favored (dilute limit), producing d — wave superfluidity 
for U <Uc^ 4.58t. 

In the optical lattice, it should be noted that the critical temperature to reach 
the d-wave superfluid state is very small since <S A;, ~ O.OAt for [/ ~ 2i [62]. 

However, adiabatic cooling may however come to the rescue when the lattice 
is gradually turned on (the system uses constant entropy trajectories). An 
estimate of the entropy of the ordered state near Tc has led to the prediction 
that one needs T/Tp ~ 0.01 to observe the d-wave (plaquette) superfluid [62]. 
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5.3 Noise Correlations and Entanglement within a BCS pair 

ARPES meaurements in condensed-matter systems allow to directly access 
the electron spectral function. This technique has played a key role in reveal- 
ing the highly unconventional nature of single-particle excitations in cuprate 
superconductors [96]. In the experimental study of ultra-cold atomic systems 
one can also resort to spectroscopic techniques, quite similar in spirit to what 
is done in condensed-matter physics. This is the case, for example, when the 
observable we want to access is a local observable such as the local density 
or the local spin density. Light (possibly polarized) directly couples to those, 
and light scattering is the tool of choice in the context of cold atomic systems. 
Recently, Stewart et al. have applied photoemission to cold atomic systems via 
an rf photon that ejects an atom from the strongly interacting system [174]. 
Additionally, Bragg scattering can be used to measure the density-density dy- 
namical correlation function [175]. In optical lattices, the lattice spacing is set 
by the wavelength of the laser, hence lasers in the same range of wavelength 
can be used to sample the momentum-dependence of various observables, with 
momentum transfers possibly spanning the full extent of the Brillouin zone. 

On the other hand, innovative noise correlations in cold atomic systems can be 
used to detect different quantum phases [176,177,178]. This technique makes 
use of the fact that the fluctuations in the momentum distribution after release 
from the trap contains information of the initial correlated quantum state. 

More specifically, what may be measured is the momentum space correlation 
function of the (fermionic) atoms in the trap: 

Gss'{k, k') = {n^pfi^,^,) - {nj{n^,p). (115) 

Such correlation techniques in expanding atom clouds have begun to be suc- 
cessfully applied in recent experiments, revealing the quantum statistics of 
bosonic or fermionic atoms optical lattices (see Ref. [55] and references therein). 

Firstly, let us imagine that atoms are released from a single macroscopic trap. 
In a typical experimental setup, the trapping potential is turned off suddenly, 
and the atoms are free to evolve independently (the atom velocity is assumed 
to be constant). This is valid if the collision cross-section is quite small; such 
conditions can be achieved experimentally. If interactions can be neglected 
during the timc-of-fiight, the average density distribution is related to the 
in-trap quantum state via [55,176]: 

{ns{f{t)))tof ~ {ns{k{r)))trap- (116) 

After a long time of flight the density distribution becomes proportional to the 
momentum distribution in the initial trapped state (n(r (t))) ~ (m/ ht) {n{k{r))); 
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the wave- vector k{'r) — mf/ (ht) defines a correspondence between position in 
the cloud and momentum in the trap. For long timc-of-fiight times, the initial 
size of the atom cloud in the trap can be neglected. In fact, in each experimen- 
tal image, a single realization of the density is observed, not the expectation 
value. On the other hand, Eq. (116) is still meaningful because each bin s in 
the image represents a substantial number of atoms Ng, while the atomic noise 
scales as 0{\/N^). In fact, since A^^, is not macroscopic the density fluctuations 
are visible. They are characterized by the correlation function: 

GsAr^r') = {nMt))nAr'm - (n,(f(t)))(n,,(P(i))). (117) 

In analogy with Eq. (116), this can be related to ground state momentum 
correlations: 

Gss' ~ Mk{r))nAk{r'))) - Mk{r))) {nAk{r'))) . (118) 

Of course, there is a proportionality constant that we ignore, for simplicity. 
Now, let us consider a superfluid state of fermionic atoms. Using the BCS 
wave-function of Eq. (75), one easily obtains that {n{r){t))tof = 2|'yg(-,-?)P [60]. 
It is interesting to observe that the latter at T = is qualitatively indis- 
tinguishable from a Fermi distribution at T = [179]. In particular, l^j^p at 
T — approaches unity well below the Fermi surface and zero well above. This 
traduces that the change in a 'metal' on cooling from TctoT — should not be 
described by the occupation numbers of one-electron momentum eigenstates. 

The essential difference between these states in fact lies in the two-particle 

— * 

correlations. For every atom with wave- vector k in the BCS ground state, there 
is another one at exactly —k. This implies prominent correlations between 
density fluctuations on diametrically opposite points in the expanding cloud: 

Gu{f,P) = 2|%>-3l>fe>^r^(r-'+r0, (119) 

where S{f + r') is a sharply peaked function of r -|- r'. This sharp peak is a 
direct analogue of the first order coherence peak in bosonic systems. It should 
be noted that the product \up]^\ is proportional to the BCS gap, ^sc- 

Now, let us consider atoms that are initially in an optical lattice. Then [176] 

{usim)) ^ (nsikim oc E e^(^-^^)-^«(4^,), (120) 

— * 

where the wave-vector k{r) now defines a correspondence between position in 
the cloud and quasi-momentum in the lattice. For example, in a superfluid 
state of bosons, where {d\dj) = \^\^ (where ^ is the superfluid order param- 

— * 

eter), n{r{t)) exhibits Bragg peaks at k{r) corresponding to reciprocal lattice 
vectors P. In the Mott state, on the other hand, {d\dj) ~ 5ij, and there is no 
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interference pattern in n{f{t)). However, it is interesting to observe that the 
D-Mott state for the plaquette system yields unusual interference peaks [180], 
which could be observed experimentally in a time-of- flight experiment. 

The density-density correlation function after a time-of-fiight takes the form: 

GsAr,r') = X; e^^-'-^>"^+^^-'-^>')(44,,d,.,,ci,,,) (121) 

ii'jj' 

+ S{f-r'){n{r{t))) - {n{r)) {n{r')) . 

In a BCS state, G|j(r, r') contains terms proportional to |A5C'P, and a super- 
fluid with \Asc\ 7^ must exhibit interference fringes at A; + A:' = nP, where 
n is an integer. The dominant contribution is for i = j and i' = j'. 

In the case of the Hubbard model with repulsive interactions, one must take 
into account the proximity to the Mott state. In the pseudogap phase, one 
expects that the phase fluctuations of the pseudogap order parameter will 
destroy the coherence peak in Gs{r, r'). In the d-wave superconducting phase, 
using the renormalized mean-field theory, one predicts that the coherence peak 
in G|j will be proportional to lAg^p; see Eq. (81). The d-wave nature of the 
state will be signaled by a modulation of the peaks with an overall enveloppe 
showing d-wave nodal directions along = -^ky and A;^ = '^k'y. In the d-wave 

BCS state of the checkerboard system, G|j,(A;, k') exhibits interference fringes 
at A; -|- A;' = nK where K is the reciprocal lattice vector of the plaquette array 
(which is half of the reciprocal lattice vector of the underlying lattice) [62] . 

Noise correlations can also be used to detect superfluidity in ladder type sys- 
tems. For example, for the spinless Hubbard ladder in the p-wave supercon- 
ducting state, using Eq. (30), we predict that G(g, q') o:\q + there is 
a singularity aX, q' — —q since 7 > 1 even though phase fluctuations cannot be 
completely avoided in ID. For electrons with spin, the singularity is more pro- 
nounced G^i oc \q + q'\^~^ and ^1/2 (Eq. (11)). For the ID case, q = k — kp 
and q' = k' + kp denote wave- vectors (momenta) relative to the respective 
Fermi points. For a comparison, it should be noted that for a single Hubbard 
chain with attractive interactions [181] C-fj, (?,?') oc \q + q'\^^^''~'^, where Kc 
is the Luttinger parameter in the charge sector (the system exhibits only a 
divergence for attractive interactions K^. > 1). In contrast, a charge density 
wave or a Mott state would always exhibit an anticorrelation in at q' = q. 

It is useful to note that the noise correlations in Eq. (118) allow to confirm that 
within the BCS state "entanglement" (singlet formation) only occurs between 
each Bell pair {k, s) and {—k, —s). The amount of entanglement entropy is: 

S = -^k^^^k - (1 - ^k) ln(l - (122) 
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where zj^ = \u^\'^ = 1 — in the BCS theory {{vj^l^ represents the probabihty 
of a pair {k, ]; —k. [) being occupied). In a Fermi hquid phase, since = 
0, then we recover that there is no entanglement entropy within a pair of 
electrons {k, s) and (— fc, —s). It is relevant to note the correspondence between 
entanglement and the superconducting energy gap via the product uj^^- 

5.4 BCS-BEC crossover and pseudogap 

At this point, it is suggestive to discuss the "BCS-BEC" crossover taking place 
in harmonic traps, and address possible connections with high- Tc cuprates 
concerning the existence of a pseudogap phase. A discussion on the BCS- 
BEC crossover can also be found in Ref. [55]. What makes those gases so 
important is their remarkable tunability and controllability. Using a Feshbach 
resonance, one can tune the attractive two-body interaction from weak to 
strong, and thereby make a smooth crossover from a BCS superfluid to a 
Bose-Einstein condensation (BEC) [182,183]. The field of BCS-BEC crossover 
is built around early observations by Eagles [183] that the BCS ground state is 
much more general than was originally thought. If one increases the strength of 
the attraction and self-consistently solves for the fermionic chemical potential 
fj,, this wave-function will correspond to a more BEC-like form of superfluidity. 
In the weak-coupling regime, fi = Ep where Bp is the Fermi energy and BCS 
theory results. However at sufficiently strong coupling, begins to decrease, 
crossing zero and then ultimately becoming negative in the BEC regime. For 
negative /x, the Fermi surface is gone and the material is bosonic. Knowing 
the ground state what is the nature of superfluidity at flnite T? Nozieres and 
Schmitt-Rink were the flrst to address non-zero temperature [184] . 

In fact, one can make relevant observations [185]. 

As we go from BCS to BEC. pairs will form above without condensation. 
In the normal state, the system takes advantage of the pairing attractive 
interaction. Only in the extreme BCS limit do pairs form exactly at T^. 

The fundamental entities in these superfluids are fermions. Pairs of fermions 
form a boson. Similar to the pseudogap phase of high-Tc cuprates, one mea- 
sures these bosonic pairs indirectly via the fermionic gap parameter A(T). 

There will be two types of excitations in these BCS-BEC crossover systems. 
Only in strict BCS theory are the excitations of solely fermionic character, and 
only in the strict BEC limit will the excitations be solely of bosonic character. 
In the intermediate ("unitary") regime the excitations consist of a mix of both 
fermions and bosons below a crossover temperature denoted T* (by analogy 
to the cuprates) while the superconducting transition takes place at Tc <T*. 
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Fig. 18. Sketch of excitations in the BCS, unitary, and BEC regimes. The yellow 
disks represent fermionic excitations. Pair excitations become progressively domi- 
nant as the system evolves from the BCS to BEC regime. 

An approach to address the BCS-BEC crossover is through the BCS Leggett 
T-matrix theory [185]. One considers coupled equations between the particles 
(with dressed propagator G) and the pairs with propagator t{P). The pairs 
and the fermions are coupled: at the mean-field level, the attractive BCS 
interaction term U between fermions produces a coupling between the pairs 
and the fermions. At the mean-field level, a pair can be converted into two 
fermions (Andreev type scattering). Pair-pair interactions are only treated in 
a mean-field averaging procedure. Below Tc, then there are two contributions 
to the full T-matrix, from the non-condensed pairs with non-zero momenta 
or pseudogap (pg) and from the condensate (sc), t = tpg + tgc where tsdQ) = 
— A|(^(5(Q)/T and we denote Q = {q,iujn) and write J2q = TY^i^^ where 
uin is the Matsubara frequency. The electron self-energy then reads: 



T.{K) = S,,(i^) + llpgiK) =Y.t{P)G,{P - K), 



(123) 



where Go is the free electron propagator: Gq^ = iun — and ~ ~ 
where eg = h^k"^ / {2m). The label pg serves to label the "pseudogap" and the 
corresponding non condensed pair propagator. It follows: 



^SC 



OJ + ei 



(124) 



This is the usual BCS self-energy. Now, one needs to include the effect of the 
non-condensed pairs. An approximation below consists to replace [185]: 



Sp,(ir) ^ -G^{-K)A. 



2 

pg- 



(125) 



This is valid below since there is a vanishing chemical potential for the 
pairs, which means that tpg{Q) diverges at Q = when T < T^. It is reasonable 
to assume that noncondensed bosons in equilibrium with a condensate must 



necessarily have a zero chemical potential. One gets A 
(this is the pseudogap equation), and the pair propagator obeys 



P9 



u 

1 + Ux{P) 



(126) 
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where x(-P) is the pair susceptibihty defined as [185]: 

x{P)-l{XuiP) + Xl^{P))- (127) 

Levin et al. assume that the latter is given by the product of one dressed and 
one bare Green's function, 

Xn(P) = E^oT(i^-m(X). (128) 

K 

(There are more than one possible choices in the normal phase.) Within this 
definition, one obtains consistent answers between T-matrix based approaches 
and the BCS-Leggett ground state equations. We have the usual BCS-like form 
for the electron self-energy ^{K) = ^j^zirji below with 

A'{T)^Al{T) + Al{T). (129) 

Since these finite momentum pairs do not have phase coherence, their con- 
tributions lead to a pseudogap in the fermion spectrum. In this context, the 
gaps add in quadrature as a result of the (valid) approximation in Eq. (125). 

It should be noted that this form of self- energy is distinct from that derived 
by Yang et al. [48] for the underdoped cuprates in Eq. (97); see Sec. 4.6. 

Knowing the electron self-energy, one finds G and tpg. In fact. Levin et al. 
argue that Ap^(T) = Z^^ J^qyto^i^g^T), where is a residue, b{x) is the 
Bose-distribution, and Qq = (f /2M* is the pair dispersion where M* is the 
"effective" pair mass [186]. It should be noted that within this approach, the 
collective modes do not couple to the pair excitations: this leads to a (f /2M* 
form of the pair dispersion. In this sense, all inter-boson effects are treated in 
a mean-field sense and (only) enter to renormalize the effective pair mass M*. 
The role of these modes at unitarity is currently unresolved. 

At T = 0, Apg vanishes since all pairs condense. 

The BEG condition Upair = gives the usual BGS gap equation A ~ A^c: 

1 + ux{^) =0 = 1 + uY.—^^^i^—- (130) 

The number of fermions determines the chemical potential n = 2J2kG{K). 
Above Tc, one must set Age — and Upair 7^ 0. By introducing the Bose- 
Einstein distribution function b{x) for non-condensed pairs, one predicts 

Alg^Z-'Y.K^1-I^Pair,T), (131) 
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Moreover, in the normal state, one gets t'g — + x(0) = Zjipair- In the 
pseudogap phase, then one can determine ^pair^ ^pg^ ^"^^ A* t)y adding the 
density equation n = 2 '}2k G{K) [187]. It should be noted that the pseudogap 
physics in cold atomic systems has been observed, e.g., in Ref. [188]. 

It is certainly interesting to note the intersection of two different fields: high- 
Tc superconductivity and superfluidity in atomic traps. In both cases, the 
pseudogap in the normal state consists of noncondensed fermion pairs. On 
the other hand, in high-Tc cuprates the pseudogap is also related to Mott 
physics which reflects the anomalously small zero temperature superfluidity; 
this aspect of the problem is not present in cold atomic systems with attractive 
interactions. Moreover, it is relevant to underline that for one particle per site 
on a bipartite lattice, the Hubbard model with U > Q maps onto the same 
model with [/ < under a particle-hole transformation; therefore, the spin 
density wave (weak couphng) regime corresponds to the BCS one and the 
Heisenberg (strong coupling) regime to the BEG one. 



6 Conclusion and Discussion 

First, we emphasize that superconductivity close to the Mott state appears 
in a variety of single-band Hubbard systems with repulsive interactions such 
as ladders, the checkerboard Hubbard model, and the homogeneous Hubbard 
model in two dimensions. The antiferromagnetic fluctuations at a large wave- 
vector (vr, vr) unambiguously favor a d-wave type singlet pairing (see Eq. (2)). 

It is maybe important to underline that the concept of an electron-electron in- 
teraction mediated by magnetic spin fluctuations is not speciflc to the cuprates 
[63] . A similar process is believed to help facilitate p-wave spin-triplet pairing 
{L — 1, S — 1) in superfluid ^He, and to lead to other pairing states in organic 
superconductors as well as heavy fermion systems [189,190]. 

Now, we summarize the main ideas of this review and address a few questions. 

Theoretical Proof of d-wave superconductivity: For quasi-lD systems, one can 
rigorously show that the (4-band) antiferromagnetic processes of Fig. 9 are 
responsible for the Cooper instability close to half-fllling (even though phase 
fluctuations cannot be avoided even at zero temperature due to the reduced 
dimensionality). For the 2D Hubbard model in the strong coupling limit, i.e., 
for the t — J model, at a mean-fleld level, one can always decouple the Heisen- 
berg exchange interaction as a BCS pairing channel, as shown in Sec. 4.5. The 
renormalized mean-fleld theory plus variational Monte-Carlo methods applied 
to the t — J model provide a good understanding of the d-wave BCS state 
for < 5 < 5c where 5c ~ 0.3 and J/t 0.2; results in the superconducting 
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state are in agreement with experiments. It should be emphasized that a single 
square plaquette (the minimal 2D unit) already sustains d-wave pairing. 

Introducing Gutzwiller-type projected BCS wavefunctions: The discovery of 
high-Tc superconductivity has certainly shown the way of using Gutzwiller- 
type projected BCS wavefunctions as trial wavefunctions to study complex 
strongly correlated electron systems. In general, the projected BCS wavefunc- 
tion is appropriate to describe a resonating valence bond insulator, a Fermi liq- 
uid phase and the superconducting region. Strong correlation close to the Mott 
state is the driving force behind the phase diagram of high-T^ cuprates and 
the resulting BCS state can be captured through a projected BCS wavefunc- 
tion. Due to the proximity of the Mott state, the situation is well distinguish- 
able from the conventional BCS superconductors which are controlled by the 
electron-phonon coupling [191]. In particular, using the standard Gutzwiller 
approximation, one gets that all the quantities containing a single-particle op- 
erator are renormalized with a factor gt vanishing as 6 at small doping levels. 

Tc.' Concerning the typical value of the critical temperature Tc in the 2D Hub- 
bard model for weak interactions, since the Cooper interaction is a 'marginal' 
interaction in the RG sense, this results in a that varies as exp —vp/U, and 
therefore is very small. For very large interactions, the renormalized mean- 
field theory predicts that close to optimal doping Tc ~ gtJf^ where J is 
the Heisenberg superexchange (we have used the fact that the variational 
d-wave gap decreases linearly with 6). A typical value for the cuprates is 
J~O.13eV0; therefore, around optimal doping where 5 ~ 0.15 and gt ~ 0.26, 
gives Tc ~ 169K. Dynamical Cluster approaches provide the same typical val- 
ues for Tc [54,192]. Cuprates close to optimal doping are governed by a Tc 
from 90K (for YBCO) to 138K (for a Hg copper oxide superconductor). The 
connection between the spin gap and J has been discussed experimentally, 
e.g., in Ref. [193]. At low doping levels, following the argument by Lee and 
Wen [102] for the destruction of the d-wave superconductivity by the quasi- 
particle excitations, one gets that Tc follows the superfluid density at T = 0, 
Ps{T = 0). From the renormalized mean- field, then one gets ps oc 5 at low 
doping levels, reproducing the Uemura et al. plot, which draws a precise rela- 
tionship between Tc and carrier density. Recently, it has been suggested that 
significant enhancement of Tc is possible in cuprate heterostructures [194] un- 
der realistic conditions providing a possible explanation for measurements on 
LSCO bilayers [195]. 

Single-Band Model: At this point, it is relevant to discuss the validity of the 
single-band Hubbard to describe cuprates. At stoichiometry the S" = 1/2 spins 
of the Cu^"*" ions in the Mott insulator order classically below a Neel temper- 



™ The relatively large Tc in the cuprates certainly results from this large value of 
J; for a comparison, note that the Debye temperature for Cu is only ~ 315K. 
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ature ~ 300K; T]v is small as a result of the small interlayer coupling. The 
highest occupied orbital is is the d,j.2_y2 orbital. On the other hand, the copper 
is surrounded by six oxygens in a octahedral environment. The Cu orbital is 
singly occupied while the p orbitals are doubly occupied. The d and the p 
orbital form a strong covalent bond resulting in a large hopping integral tpd- 
Using the hole representation, the oxygen is empty of holes and lies at an en- 
ergy Ep which is higher than Ed- In the parent compound, the virtual hopping 
process tpd then induces an effective exchange interaction: J = t^^/ {Ep — Ed)^ 
where the effective hopping amplitude t ~ t^J {Ep — Ed) and Ep — Ed plays 
the role of the Hubbard U in the one-band model. Experimentally, an energy 
gap of 2eV is observed from optical conductivity. By doping the Cu-0 planes, 
a doped hole resonates on the four oxygens surrounding a Cu and the spin 
of the doped hole combines with the spin on the Cu to form a Zhang-Rice 
singlet [196]; the Zhang-Rice singlet can hop from site to site with an hopping 
amplitude ~ justifying the single-band Hubbard model. This effectively pro- 
duces sites with a Cu^"*" valence. A tight-binding model for Cu can be built 
thoroughly; see Ref. [34,90]. 

Pseudogap phase: Even though the d-wave superconducting state of the cuprates 
can be relatively well-described through a projected BCS wavefunction, the 
pseudogap phase still remains a "delicate" issue from the theory side. Several 
theories predict the possibility of an RVB-type state, i.e.., the existence of a 
pseudogap phase close to half-filling with preformed pairs and a truncated 
Fermi surface. More precisely, the ladders are nice prototype systems allowing 
to show the possibility of a truncated Fermi surface close to the Mott state 
[30] where antinodal regions form the celebrated D-Mott state, which is char- 
acterized by a spin and a charge gap (see Sec. 3 and 4.2). The ladder systems 
also provide a quantitative expression for the electron's Green function in the 
pseudogap phase and allow to prove the (generalized) Luttinger Sum Rule. 

In two dimensions, the RG approach in the 2-patch model, when the Fermi 
surface lies at the saddle points (0, tt) and (tt, 0) (for t' /t ~ 1/4 this happens 
for 5 ~ 0.2), allows to predict a novel phase close to half-filling in the t — t' — U 
model in which the antinodal directions become incompressible as a result of 
elastic umklapp processes. The opening of a spin gap has also been confirmed 
through numerical calculations at the low-energy fixed point [46] . This incom- 
pressible spin-gapped region yields some resemblance with the D-Mott spin 
liquid region of the ladder system. By increasing the electron density, from the 
quasi- ID analysis for t±_ ~ t, one predicts the propagation of this spin-gapped 
Mott region along the umklapp surface; this weak-coupling RG approach also 
suggests the appearance of Fermi arcs in the pseudogap phase. In the strong- 
coupling limit, the SU(2) slave-boson approach developed by Wen and Lee [51] 
also tends to predict the existence of Fermi arcs along the nodal directions, 
through a spinon-holon recombination. However, several approximations are 
used to study the effect of gauge fiuctuations on the electron Green's function. 
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Finally, by combining the results from the ladder systems [30] and those from 
the renormalized mean- field theory [12], Yang et al. [48,49] have proposed an 
ansatz for the single-particle Green's function in the pseudogap and in the 
superconducting phase, which shows very good agreement with ARPES and 
STM measurements. A comparison with other existing theories for the electron 
propagator in the pseudogap phase can be found in Ref . [137] . 

Nodal- Antinodal Dichotomy and 2 gaps versus 1 gap: It is relevant to observe 
that the RG approach in the weak-coupling limit and the slave-boson approach 
in the strong-coupling limit suggest two distinct scenarios for the occurrence 
of the pseudogap in the antinodal regions. The RG approach in the t — t' — U 
model shows the occurrence of an incompressible phase for U > Uc where 
Uc/t oc ln~^(i/f'). As a reminiscence of the ladder systems, this pseudogap 
phase then can be seen as a precursor of the Mott transition. Below T^, for 
the underdoped regime, superconductivity proliferates in the Fermi arc regions 
only as a result of the antiferromagnetic fluctuations (preformed pairs) in the 
antinodal regions; this scenario suggests that superconductivity in the nodal 
regions appear through an Andreev type process involving the hot (antinodal) 
regions. The quasi- ID and the weak-coupling RG approach in two dimensions 
predict the appearance of two energy scales and two distinct gaps close to 
half-filling: T* embodies the energy scale at which umklapp scatterings play 
an important role in the antinodal regions opening a charge and spin gap. Tc 
represents the superconducting energy scale at which the Fermi arcs become 
superconducting. One important consequence is that the energy gap in the 
nodal regions does not scale with the pseudogap at (tt, 0). The two gap scenario 
is supported experimentally, by ARPES measurements, Andreev reflection 
studies, and Raman scattering. This is contrast with the Bose condensation 
scenario which only involves a single gap, i.e., the d-wave pseudogap. A direct 
evidence for competition between the pseudogap and superconductivity has 
been recently reported using ARPES [197]. 

Fermi liquid for severely overdoped systems: The Fermi liquid phase is rela- 
tively simple to capture. From RG arguments, it appears when antiferromag- 
netic fluctuations and umklapp scatterings are negligible due to the relatively 
large doping. The Cooper processes then become marginally irrelevant. In the 
large U limit, the antiferromagnetic fluctuations are governed by the Heisen- 
berg exchange coupling J. The emergence of the Fermi liquid phase can be 
understood from the competition between the exchange energy J and the ki- 
netic energy which is of order t per hole or 5t per unit area. When St ^ J, 
we expect the kinetic energy to win and we expect a Fermi-liquid metal with 
a weak residual antiferromagnetic correlation. This fact is also transparent in 
the slave-boson theory where spin and charge bind when the d-wave gap is 
zero. One expects that the crossover to the Fermi liquid phase should occur 
when (5i ~ J or 5 ~ J/i ~ 1/3 (in agreement with experiment at T — > 0). 
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Breakdown of Landau-Fermi liquid in the overdoped regime: When St < J, in 
general, the outcome is less clear because the system would like to maintain 
the antiferromagnetic correlation while allowing the holes to move as freely 
as possible. In particular, the evolution of (or to) the antiferromagnetic state 
at low dopings 5 <S 1/3 remains a subtle question due to the proliferation 
of alternative states (d-density wave state, stripe phases, for example) even 
though the RVB theory emerging from the ladder analysis gives a reasonable 
(and quantitative) idea of the pseudogap phase. It is interesting to observe 
that Abdel-Jawad et al. [138] have reported a flagrant correlation between 
charge transport and superconductivity in overdoped high-Tc cuprates. Re- 
cently, Ossadnik et al. [149] have shown that d-wave pairing in the overdoped 
region of the phase diagram is strongly related to the appearance of a strongly 
anisotropic scattering vertex both in the particle-particle and particle-hole 
channels. Through a functional RG analysis, they have found a scattering 
rate with an isotropic Fermi-liquid like part and a strongly anomalous part 
that varies linearly with T and presents a d-wavc character. This scattering 
rate with d-wave symmetry emerges as a result of antiferromagnetic processes 
connecting the antinodal points and must be attributed to the increase of the 
4-point vertex when decreasing the energy. This phenomenon can be perceived 
as a precursory effect of the pseudogap region (which should emerge as soon as 
the Fermi surface touches the saddles points producing a van Hove singularity 
and an incompressible spin-gapped antinodal region as soon as U > Uc-) 

Other ingredients? There are many reasons to believe that the single-band 
model contains most (but maybe not all) of the ingredients necessary for un- 
derstanding high-temperature superconductivity in the cuprates. This does 
not preclude that additional phenomena, such as interactions between elec- 
trons and the lattice (phonons) and inhomogcncitics at the mcsoscopic level, 
must also be included to have a quantitative description of the phase di- 
agram. In particular, recently it has been proven that the stripe order in 
La2_5Ba5Cu04 frustrates the three-dimensional phase order, but is fully com- 
patible with 2D superconductivity [119, 120]. The strict interplanar decoupling 
arises because the planar superconductivity contains a periodic array of lines 
of TT-phase shift which rotate of 7r/2 up the c-axis together with the spin and 
charge stripe ordering in the low-temperature tetragonal phase. 

Why is the high-Tc problem hard, hut important? The superconducting region 
close to the Mott state may already be unconventional and is an interesting 
question on its own. The weak-coupling RG approach in the quasi-lD limit 
and for the t — t' — U model in two dimensions suggests a d-wave supercon- 
ducting phase (with two distinct gaps) where only the Fermi arcs participate 
in the superfluid weight, whereas the antinodal regions remain incompressible. 
By increasing the doping further, the superconducting region then proliferates 
along the antinodal directions. The coexistence between Mott and supercon- 
ducting regions in momentum space is an interesting idea that has to be pushed 
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further. The high-T^; problem is hard because one must take into account the 
proximity to the Mott state that produces prominent phase fluctuations: the 
Gutzwiller projection is difficult to implement beyond the renormalize mean- 
field theory and the slave-boson method suffers from several approximations in 
the treatment of gauge fluctuations. Understanding the anomalous pseudogap 
phase in the underdoped cuprates remains a central challenge. Rather than 
deal directly with the very difficult strong correlation, one can relax the condi- 
tion and approach the problem assuming a weak to moderately strong onsite 
interaction. In fact, in the quasi-lD regime and in the 2-patch model, the RG 
approach gives quantitative (and encouraging) results close to the Mott state. 
Prom an experimental point of view, the situation is also encouraging: it now 
appears possible to probe the normal phase(s) at zero temperature by applying 
a large magnetic field [130]; this represents a powerful new tool to study the 
normal regime and especially the underlying physics in the antinodal regions. 

Similitudes with Fe-pnictide superconductors ? The recent discovery of super- 
conductivity in a family of Fe-based oxypnictides with quite large transition 
temperatures has led to tremendous activity aimed at identifying the pairing 
mechanism in these materials. The high transition temperatures and the elec- 
tronic structure of the Fe-pnictide superconductors suggest that the pairing 
interaction is of electronic origin. Additionally, the involvment of antiferro- 
magnetic correlations in the superconductivity also seems to emerge. On the 
other hand, a few distinctions already appear between the two classes of uncon- 
ventional superconductors: first, the Fe-based oxypnictides are described by a 
complex multi-band system and second, in the parent iron based compounds, 
antiferromagnetism seems to emerge from a collective spin-density-wave order 
instability of an itinerant system (and not from local moment physics) . 

Ultracold atoms as Quantum Simulators: Ultracold atoms in an optical lattice 
allow in principle to engineer the Hubbard model in a controlled way and 
to detect d-wave superfluidity via novel tools such as noise correlation mea- 
surements [60]. In this sense, this would represent a quantum simulator as 
Feynman first envisaged it. Spectroscopic tools have been successfully imple- 
mented and it is now possible to apply photoemission or Bragg scattering type 
experiments to study cold atomic fermions. In fact, there are many promising 
roads for research with ultra-cold fermionic atoms. In particular, optical lat- 
tices have opened the possibility of studying strongly correlated regimes, such 
as high-temperature superfluidity and its interplay with Mott localization. 
They also offer the possibility of studying these systems in regimes which are 
not usually reachable in condensed-matter physics {e.g., under time-dependent 
perturbations bringing the system out of equilibrium), and to do this within 
a highly controllable and clean setting. At a very general level, optical lattices 
allow to engineer the many-body wave-function of large quantum systems by 
manipulating atoms individually or globally. On the other hand, to probe the 
Neel ground state of the half-flUed Hubbard model and the d-wave superfluid 
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state when doping, this requires to reach initial temperatures in the trap such 
that the entropy per particle is smaller than In 2 (we have set ks = 1)- It would 
also be worthwhile to understand more deeply the physics of atomic fermions 
at the unitary limit and the pseudogap phase in these systems. Finally, it 
is also relevant to underline that optical lattices also offer the possibihty of 
exploring exotic forms of superfluidity out of equilibrium [198]. 

Acknowledgments: We thank our collaborators E. Dagotto, N. Furukawa, 
C. Gros, C. Honerkamp, J. Hopkinson, R. Konik, A. Lauchli, U. Ledermann, 
M. Ossadnik, M. Salmhofer, M. Sigrist, A. Tsvehk, M. Troyer, K. Y. Yang, 
and F. C. Zhang. K. L. H. is also grateful to the National Science Foundation 
(NSF) for financial support related to this research, under the contract DMR- 
0803200, and to the members of CIFAR for stimulating discussions. 



A Renormalization Group 

Renormalization Group (RG) techniques have been developed in particle physics 
to investigate the high-energy behavior. In condensed matter systems, RG 
techniques are rather applied to study the low-energy condensed matter physics. 

The Kadanoff- Wilson type high-energy mode ehmination (see e.g. the review 

by Shankar [87]) has been successfully applied to describe the physics of a 
single Hubbard band. The kinetic Hamiltonian Hq describes Dirac fermions 
with the velocity v = 2tsin(A;i?) and the interaction part Hjnt takes the form: 

/ dk \ 

Hint =[Ilj_^ 2^r^^' + k2-k3- k^)hint, (A.l) 

where 

hint^9i'^Rsiki)'^{s'ik2)'^Rs'ih)'^Ls{ki) 

+ (^l) ^L' (^■■2) *L.' {h)'^Rs{h) 

where we have introduced s =J, if s =] and vice- versa. Here, we have neglected 
the umklapp scattering assuming that the system is not half-flUed. The g4 
term includes completely chiral interactions. This term does not affect the low- 
energy properties, but only renormalizes the Fermi velocity (this is transparent 
when using bosonization); therefore in this work we systematically drop such 
terms. The 512 term corresponds to forward scattering whereas gi describes 
backward scattering (or Cooper scattering). 



(A.2) 
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The Kadanoff- Wilson scheme is as follows. Let us denote S the action corre- 
sponding to H = Hq + Hint and to the cutoff A, and the action corresponding 
to the kinetic term (Dirac fcrmions) reads: 

p=HI L,s 

where ± corresponds to right and left movers respectively. The action S' of 
the "energy-reduced" system is given by 



" = j D^>e-^, (A.4) 



where the (Grassmann) integration is carried out over fields with A/6 < 
|A;| < A and 6 > 1; the wave-vector k is measured from kp or —kp- The 
renormalization group transformation is completed by a rescaling procedure 
which returns the cutoff to its original value (r = it is the imaginary time): 



X 



bx;r ^br;^ ^b-^/^^, (A.5) 



such that the action S'q — 5*0 (A/6) = 5*0 (A) is invariant under the RG proce- 
dure. Now, let us decompose exphcitly ^ — ^ + ^', where ^' is the field to 
keep: 

where 5"/^ denotes the interacting part containing j times the field ^. More 
precisely, we have decomposed the noninter acting action as 5*0 (^P) = 5*0 (^') -|- 
S'o(^) = Sq + Sq. The interacting part is a bit more involved and one has: 

Sj = Sio + Y.Sij. (A.7) 

The integration over ^ is not straightforward and one has to perform a pertur- 
bative expansion. The important contribution stems from the second order: 

S'^S', + Sio-l/2{S]2)-o, (A.8) 

where (...)o indicates averaging over ^. The calculation is done diagramati- 
cally. Taking into account particle-particle and particle-hole diagrams (that 
differ from an overall minus sign), in the action S', the couplings gi and g2 
are renormalized as: 

r 2 1 r In 6 2 / a o\ 

dgi^ c/i and = -7^ — gi- (A.9) 

TTV ZTTV 

Iterating the RG process A — > A/6, A/6 — > A/6^,..., then results in the RG 
equations for the couplings, 

= and = 1 (A.10) 

dl TTV dl ZTTV 
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where I — \nb. This shows that gi — 2g2 — const, and that 



gJl) ^ 3]^^} (A.ll) 

For repulsive interactions, 5'i(0) = (72(0) = f/ > 0, then gi vanishes in the 
low-energy limit and g2 = U/2 in agreement with the Luttinger liquid. In the 
text, for a given band, say 1, we use the notations cj^ = —{gi/2 — g2) > and 
Cn = 2^1. 

For the N-lcg ladder case, the Operator Product Expansion (OPE) allows to 
derive the one-loop RG equations in a quite simple way. The RG equations 
for the situation close to and at half-fiUing are given in our Ref. [30]. This 
approach systematically takes into account the U(l) symmetry for charge and 
SU(2) symmetry for spin. The charge and spin currents have been introduced 
in Sec. 2.3. A product of normal-ordered currents obeys the following short- 
distance expansion for z ^ (Wick-Theorem): 



jRiz)jR{0) = : ^i,{z)^ns{z) :: *L'(0)*i?.'(0) : (A.12) 



27rz ' ^ ' 2t:z 

here, we have denoted z = x + ivr for right movers (for left movers z — > 
z*). Now, when including explicitly the different bands, one gets the (non- 
umklapp) currents (for convenience, we redefine: Zi — ViT — ix): 



jRij (A. 13) 



jRijJmm ~ 


Sjl 

2ttzj 


J Rim 


2'KZi 


ja jb 
'^Rij'^Rlm ^ 


^ab , 

T I 


27r Zj 


J Rim 


jRijjRlm ~ 


2ttzj 


ja 


27rzi 



ja 

where we have introduced the totally antisymmetric tensor e"^'^. These rules 
can be easily extended by including umklapp scatterings. Since all the cou- 
plings are marginal, the RG equations take the general form: 

dgi 

— = Aijkgjgk- (A. 14) 

For the (spinless) two-leg ladder, it is straightforward to show that the inter- 
band Cooper interaction C12 will renormalize all the other coupling channels. 
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It is also instructive to fix the Hubbard initial values of the various couplings; 
to avoid any mistake, one can make use of certain rules, such as: 

jRijJum - ^JRij ■ Jum = 2(^'L^Ri.^L-^Lm.- " ^'L^Ri^^Ls^Lms) • (A.15) 

In the case of spinful fermions, this results in: 

3U 

"^fkh ^ = 2cIj^ = _^ -j^-j ) (A. 16) 

where {k — l,k — N), [k — 2,k — N — 1),... , and for i^k,k one gets: 

fik — ^ik — 2(]\r+T) ■ (A- 17) 

For the SU(2) spin interactions one gets /4 = ^fik and c-^ = 4c^^ for all i and 
k. To determine the bare values of the umklapp processes, one can resort to: 

liijlum = -^ijs^Lls'^^ms^Lm-s + "^^s^ Lms^^Rii^ Us. (A.18) 

One also infers that: 

'^kkkk ^ ^'^kkkk ^ 2(_/V + 1) ' (A. 19) 

whereas w^^^^^. = at the bare level. 



B 4-band interactions and quasi- ID antiferromagnetism 

B. 1 Definitions 

Here, let us assume that the number of legs is even. 

The 4-band interactions combine processes between two different band pairs 
{k,k) and {1,1) (Fig. 9): 



N/2-1 N-l 

E E 

1=1 k=i+i 

where, 
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N/2-1 N-l 

Hab^ E y^^C (B.l) 



^Ik = Ci^kTiJlilkJLkl + JRlkJLkl + h.c.) (B.2) 

~ ^ikki {jRik • Jluc + Jm ■ Jlh + h.c)j 

+ <kM (iRikhu + Im^Lik + h.c) 
+ ^zkfc ijmihkk + ^RkkhiT + ^-c-) 

Since t;;;; = this also implies: 



C,S C,S /-p) Q\ 

^lkkT~^lkkT \^-'->) 

The bare values of the four-band interactions: 

The 4-band couplings describe antiferro magnetic processes, i.e., terms such as 
'^Irs'^ils' involves a longitudinal momentum transfer of tt. 

The RG equations are given explicitly in our Ref. [30]. The key point is that, 
for a small number of chains, say N — A, the 4-band interactions stay in a 
fixed ratio [30]: 



^1234 ' U 



5/12 



oc - -0. (B.5) 

•-14 X^U/ 

At low energy, the bands 1 and 4 are therefore decoupled from the bands 2 and 
3. The low-energy Hamiltonian is thus the sum of two-leg ladder Hamiltonians 
H — i?i4 -|- H23. All charge and spin excitations are therefore gapped (similar 
to the Heisenberg 4- leg spin ladder) . 

For N odd, there is an umklapp term for the band / = {N + l)/2, 

I dxu'ii(lRiilLii + h.c.), (B.6) 
leading to a Mott gap in this band as well. 



B.2 Quasi- ID antiferromagnetism 



Now, we study the quasi-lD antiferromagnetic phase. The RG equations for 
the 4-band couplings in the (simplified) limit of large N are given in Ref. [36]. 
We have checked that the 4-band couplings converge to the precise values [36] 

t ~ Sgik = 4c'{,-,i = ^c^kkT = ^^IkkT = -^<-kT- (B.r) 
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Additionally, one can check that wj^-^j^ remains small. For large A^, the 2-band 
and single-band interactions are dominantly renormalized by the 4-band in- 
teractions. The asymptotic ratios of the 2-band and single-band interactions 
can be calculated by inserting the ratios (B.7) in the RG equations for these 
couphngs [30]. One then finds that the couplings of the band pairs {k, k) fiow 
and approach (novel) fixed ratios [36]: 

t - Sg,, = Afl-, = 3/4 = 4c4 = 3cU = Aul^,, = 8ul,-,-, = 3«4^,. (B.8) 

The main difference with the D-Mott phase is that now flows to strong 
coupling and that clf^ remain small (but attractive). For bands k and / that 
are close together on the Fermi surface, k ^ I, we check that the 4-band 
couphngs become the same as the corresponding 2-band couphngs. In the limit 
— > oo, the RG equations simplify considerably and an analytic solution is 
in fact possible. In particular, the energy scale at which the system enters in 
the antiferromagnetic phase and opens a uniform Mott gap can be determined 
using that [36]: 

d 1 

■j^sik = —siiSki, (B.9) 



where 



si^='^-^ + c^,M + <M- (B.IO) 
For the extreme case t± — > t, the scale of divergence is given by: 

where U < t/lnN for the validity of our calculations. For t±/t — 0, we find: 

^c=|. (B.12) 

The logarithmic correction for t± ^ t stems from the van Hove type singularity 
(fi = vn ~ t/N)- This allows us to check that a quasi- ID antiferromagnetic 
type ground state occurs at the energy scale 

Eafm ~ te-^' = te-"^'^, (B.13) 

where c is a function of t±/t. Note that the 4-band antiferromagnetic processes 
with a weight oc A^ are responsible for the RG instability at the energy scale 
Eafm] see Eq. (B.9). This also favors the occurrence of a uniform Mott gap. 

Remarkably, with the ratios of the couplings given above, one can check that 
the low-energy Hamiltonian for energy Ec{N) < E < Eafm takes the form 

H = Ho-^Y,J 9ijT^i ■ T^j, (B.14) 
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where Hq is the kinetic term and the staggered magnetizations obey: 



K = ^Ls<s'*fcLs' + h.c. (B.15) 

In fact, the 2-band couphngs f^^ and m^^^^ give the contributions rhk -rhk and 
the 2-band couphngs c^'| and m^'I^^ lead to the products mfc • m^. The 4-band 
couplings provide the other couplings in rhi ■ rhj. For large A^, one can check 
that all couplings take the same value gij = g > 0. This represents an effective 
spin Hamiltonian for the half-filled weakly-interacting Hubbard Hamiltonian. 

As a result, all the band pairs {k, k) are interacting with each other. In par- 
ticular, for N odd, the (nodal) band j = {N -\- 1)/2 is now interacting with all 
the other band pairs and there is no fundamental difference between odd and 
even N. This contrasts the ladder-case at energies smaller than Ec{N), where 



only interactions within the band pairs (fc, k) are present] I: this leads to an 



odd-even effect in the spin sector that is reminiscent of the spin ladder. 

One can also resort to bosonization to analyze the low-energy properties [36]. 
Interestingly, the charge sector is reminiscent of the two- leg ladder system, i.e., 
for a band pair {k, k), the total charge mode 4>kkc+ is pinned to zero whereas 
in the antisymmetric charge sector 9p.j^^_ ~ 0; see Sec. 2.3. In contrast, in 
the spin sector, both the fields 4>kks± and the dual fields Oj^i^^ appear in a 
cosine, resulting in a competition between different phases. However, since 
the Luttinger parameters in the spin sectors (for a given band pair [k, k)) now 
obey > 1 and Kj^ig_ < 1, this favors the pinning of the fields (p^kg- 

Oj^i^g_^_. This is the main difference with the D-Mott phase where one rather 
gets Kf^igj^ < 1 and Kf.j.g_ < 1. 

To be more precise, only the differences 0^^^+ — Ongj^ appear in the 4-band 
Hamiltonian and as a result the total magnon mode, given by 

= y ^ 51 ^k-ks+ and 0T = y ^ 5^ , (B.16) 

remains gapless [36]. The real-space spin operator at the site x in chain i reads: 

^[(^) = \lll^klfk (^Uo-L'^SLs'e^^'^'^^'-^^^ + R^l), (B.17) 
^ k 



where 7j>„ = y'2/(A^ + 1) sin(7rjm/(A^ + 1)) (consuh Eq. (49)) and at half- 
filling kpk + kpi = vr for the band pair {k, k). Using the expression of 5'[(x) in 



In particular, the 4-band processes become suppressed due to the occurrence of 
the spin gap in the band pair {k, k) or as a result of < 1 since c^^ changes of 

sign. Here, Kf^-^^^ denotes the Luttinger parameter associated with the symmetric 
spin mode within a band pair {k, k). 
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terms of the different bosonic fields, one gets: 



{S^{x) ■ Si{0)) oc i-iy^' cos 



(7ra;)/x^/^, 



(B.18) 



and X < l/Ef.{N) since we are in tfie quasi-lD antiferromagnetic pfiase. 

It is relevant to remember that the quasi- ID antiferromagnetic phase has the 
same charge properties as the D-Mott state. By doping, the chemical potential 
jjL couples to the total charge mode of each band pair (l)kkc+- Since all the 4- 
band interactions contain the charge field (t)kkc+) it implies that they will be 
suppressed at low very energy scales. This shows that a (quasi-) long-range 
order is difficult (impossible) to stabilize away from half- filling, i.e., in the 
presence of mobile holes. On the other hand, the effect of doping on the 2- 
band interactions is the same as when doping the D-Mott state. 

Therefore, close to half-filling, we still expect the occurrence of a truncated 
Fermi surface where the holes enter first the nodal directions and the antinodal 
directions form a D-Mott spin liquid state. 



C Slave-boson approach for the t — J model 

Here, we briefiy review the slave-boson approach to the t — J model. For more 
details, consult Refs. [14,199]. The slave-boson approach suggests that antifer- 
romagnetic fluctuations (singlet pairs) are prominent in the pseudogap phase 
below T*. This explains the reduction of the magnetic susceptibihty and the 
reduction of the c-axis conductivity. On the other hand, within this approach 
that relies on spin-charge separation, superconductivity in the underdoped 
regime occurs when the holons (holes) condense [2]. In quasi-2D systems, the 
Bose-Einstein condensation temperature is roughly proportional to the areal 
density. Above Tc, preformed pairs exist, but the hole motion is incoherent 
resulting in "gauge" fluctuations. This approach predicts a single gap scenario 
but suffers from several approximations in the pseudogap phase. 



C.l Slave-Boson scheme and Bose condensation 

The effect of the strong Coulomb interaction is represented by the fact that 
the electron operators dis and djg are the projected ones in which double 
occupancy is forbidden. This means: 



(C.l) 



s 
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A method for treating this constraint is the slave-boson approach [200,201]. 
This method was successfully developped for the Kondo problem. The electron 
operator is then decomposed as: 



where en = — e^i = 1 is the antisymmetric tensor, f-^ is the fermion operator 
whereas 6j, dj^ are the slave-boson operators. More precisely, aj refers to a 
doubly-occupied site whereas bl corresponds to the vacant state. Within this 
formulation, the constraint becomes: 

E fUs + b% + 4d, = l. (C.3) 



In the limit of large onsite interaction, one can ignore doubly-occupied sites 
and as a result the projected electron operator is written as: 



4 = fisbi, (C.4) 



with the prerequisite: 



flfis + blk = l. (C.5) 



This constraint can be implemented through a Lagrange multiplier Aj. The 
idea is to write down the electron operator as the product of boson (charge or 
holon) and fermion which carries the spin index. 

The Heisenberg term Si ■ Sj —niUj/A where rii = (1 — can be re- written in 
terms of the spin operator fig if we neglect the nearest-neighbor hole-hole inter- 
action blbibjbj (this is reasonable if one considers cases close to half-filling with 
a small number of vacant sites). One can then apply a Hubbard-Stratonovich 
transformation to decouple the exchange term both in the particle-particle 
and in the particle-hole channel. The partition function then becomes: 

J DfDpDbDb*DXDxDA exp drL^ , (C.6) 



(formally / becomes a (Grassmann) fermionic field and b a boson field) where 
the Lagrangian obeys: 
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(C.7) 



-J 



J 



where Xij represents fermion fspinon) hopping and Aij represents fermion 
(spinon) pairing I. and J ~ J l^^l . Now, it is relevant to observe that the 
Lagrangian is invariant under a local U{1) (gauge) transformation: 



/.s-e*'7.s, (C.8) 

Note that fermions and bosons are coupled to the same gauge field. The mean- 
field theory corresponds to the saddle-point solution to the functional inte- 
graOl] idL/dx*j = and dL/dA*j = 0): 

s 

Similar to the renormalized mean-field theory, at half-filling, a variety of mean- 
field approaches give identical energy and dispersion [14], as a result of the 
underlying SU(2) particle- hole symmetry: 

4 ^ 7.4 + v^f^l h - -^;4 + (c-10) 

Essentially, in the half-filled case, this approach suffers from the same problems 
as the renormalized mean-field theory of Sec. 4.5. Note that at half-filling, 
there are no bosons and the theory is purely that of fermions (similar to the 
renormalized mean- field theory). 

On the other hand, by hole doping, one expects that the best state is the 
projected d-wave superconductor, and that the theory becomes better defined. 
Essentially, the projected staggered fiux state always lies above the projected 
d-wave state (and the energy difference goes to zero when 5^0). The behavior 

They correspond to the two ways of representing the Heisenberg exchange inter- 
action. 

A discussion on the precise value of J is given in Refs. [14,199]. 

The saddle-point value of the Lagrange multiplier is Aj = A. For the half-filled case 
A = = /i and away from half-filling from half-filling one must satisfy ^^gifisfis) = 
1-6. 
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of bosons is crucial for charge dynamics. Within the mean-field theory bosons 
are free and condensed at Tbe- The bosons see the same band dispersion 
as the fermions and condense at the bottom of the band minimum at low 
temperatures. Assuming a weak three-dimensional hopping between 'layers' 
one obtains a finite Tc ~ Tbe that is proportional to the boson density S 
[202] . Below this temperature, we have electron pairing because the BCS order 
parameter {d^-^d_^^) = b'^{f^if-ki) ^ where (bi) = 6 = v^. 

This slave-boson approach supports the scenario by Anderson [2] that spin- 
singlet formation (RVB) turns into superconductivity via the Bose conden- 
sation of holons (vacant states). Kotliar and Liu [202] and Suzumura et al. 
[203] found d-wavc superconductivity in the mean-field theory. In the d-wave 
superconducting state, one gets: {bi) 7^ 0, % 7^ 0, and A 7^ 0. In the pseudo- 
gap phase, spin-singlet formation subsists but there is an incoherent charge 
motion: A and x ^-re nonzero while 6 = 0. Within this approach, A represents 
the pseudo-gap and the superconducting gap: this is a single gap scenario. 

It should be noted that the U(l) gauge theory (including Gaussian fluctuations 
around the saddle point) suffers from several limitations in the underdoped 
regime [14]. 

C.2 SU (2) formulation 

Then, an SU(2) theory for underdoped high-T^. cupratcs taking into account 
Eq. (C.IO) has been developed. At a general level, one can formulate an SU(2) 
gauge theory by introducing the fermion doublets: 

^l-ikfli) *5 = (/a-4)- (C-ii) 

Then, one can introduce an SU(2) version of the slave-boson theory [51]: 




(C.12) 



where we have introduced the (bosonic) vector {hiY = (6|j 621) ^^^1 similarly 
for the spin (spinon) part, and ■0 = ir^ip* where represents the second 
Pauli matrix. The t — J model can now be written in terms of the fermion- 
boson fields. The local SU(2) singlets must satisfy {ip\fipi + h\fhi)\phys) = 
to ensure that the Hilbert space of the t — J model is respected (operators 
ipifipi that generate SU(2) local transformations vanish within the physical 
Hilbert space); Tj are the Pauli matrices. On a given site, there are only three 
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states that satisfy the above constraint: f^\Vac), fj\Vac), and {l/^/2){b\ + 
b2flf^)\Vac) , corresponding to a spin- up a spin-down fermion, and a vacancy, 
respectively (again, \Vac) refers to the vacuum). The constraint leads to: 

E(/i/-) + - {blb2^) = 1- (C.13) 

s 

Within the SU(2) gauge theory, the density of fermions J^sifLfis) is not nec- 
essarily equal to 1 — S, where the hole doping is defined as {bifiii + &2j^2i) = 5- 
This is because a vacancy in the t — J model may be represented by an empty 
site with a hi boson or by two fermions (spinous) in a singlet state with a 
62 boson. It describes the physical idea that adding up a spin-up fermion or 
removing a spin-down fermion give the same state (after projection to the 
subspace of singly occupied fermionsj^ 



In the superconducting phase, one must satisfy (621) = such that the SU(2) 
approach agrees with the U(l) slave-boson theory. 



The mean-field Hamiltonian takes the form: 



(C.14) 



^ ti^l^ijhj + h.c.) - fiYl hi^i 

i 



2 .. ., 



The averaged constraint (ipiTipi) = has been enforced through the Lagangian 
multipliers Oq, and oq is invariant under spin rotation, and 



U.J 



In this SU(2) formulation the d-wave superconducting phase is described by 
the following mean-field Ansatz (by analogy to the U(l) slave-boson mean-field 
theory) [51]: 

Ui,i+^ = -Xr^ + Ar^ , ^7^,^+^ = -^r^ - Ar\ (C. 15) 

al^O, a'/ = 0, (61)7^0, (62) = 0, 

In the superconducting phase, the fermion and boson dispersion are given by 



25 



Annihilating a spin-up electron is produced by applying /| onto the vacancy state. 
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±Ef and ±Eh — jj, where 
and 

3J 3J 

ef = — — {cos kx + cos ky)x, Vf = — — {cos — cos ky) A. (C.17) 

Additionally, one obtains: 



E, = ^{et + al)^ + vl (C.18) 

where 

eb = —2t{cos k^ + cos ky)x, rj^ = —2t{cos k^ — cos ky) A. (C.19) 

The bosons condense (at = 0) at the bottom of the band minimum. 

Interestingly, the Fermi liquid phase is similar to the superconducting phase 
except that there is no fermion pairing (A = 0). The Fermi liquid contains bo- 
son condensation. In this case, the electron Green's function (Sd) oc {tp'^hh^ip) 
is proportional to the fermion Green's function (ip'^ip). Thus the electron spec- 
tral function contains a 6-peak in the Fermi liquid phase. 

C.3 Spinon-holon binding in the Pseudogap phase 

Within the SU{2) slave-boson approach, the "pseudogap" phase is described 



by a staggered flux liquid (sfL) phase I. where [4] 



Ui,+, = -xr' - t{-iy^+''yA, U,,i+y = -XT'" + t{-iy^+'yA, (C.20) 

4 = 0, (&1,2) = 0. 

In the sfL phase, the fermion and boson dispersions are given by ±£^/ and 
±£"6 — /i, but now Oq = 0. Since Og = we must have J^sifLfis) = 1 cind as 
a result {blbi) = (62^2) = S/2. In fact, the sfL phase and the d-wave pairing 
phase of the fermions are connected via a site-dependent SU(2) transforma- 
tion that is exphcitly given in Refs. [14,51]. Moreover, due to the U(l) gauge 
fluctuations, the antiferromagnetic spin fluctuations in the sfL phase are as 
strong as those of a nested Fermi surface despite the presence of a pseudo- 
gap. In practice, the sfL phase resembles an exotic mixed Spin Density Wave 
(SDW) + superconducting (SC) state. 



The physics is similar to non-interacting electrons on a lattice pierced by magnetic 
flux lines whose directions alternate in a checkerboard pattern. 
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Now, it is relevant to focus on the electron spectral function in the pseudogap 
phase (that corresponds to the sfL phase in the SU(2) slave-boson approach). 
At the mean-field level, the electron Green's function is given by the product 
of the fermion and boson Green's function, so the electron spectral function 
is a convolution of the boson and fermion spectral functions. 

In the pseudogap phase above Tc, a first approximation generally used in the 
literature consists to argue that the boson is almost condensed, so the boson 
spectral function contains a peak at A; = 0; the weight of this peak is ~ 5 
and the width is ~ T. This leads to the following mean-field BCS-like electron 
Green's function: 




Gel = ^\ V + + G^nc. (C.21) 



where Uj: and are the BCS coherence factors: 



The peak in the electron spectral function crosses zero energy at four points 
(±7r/2, ±7r/2). Thus, the mean-field sfL phase has four Fermi points (similar 
to the d-wave phase). 

The second approximation consists to argue that the gauge fluctuations pro- 
duce a strong attraction between bosons and fcrmions due to the fluctuations 
around the mean-field state. This attempts to bind bosons and fermions and 
form Fermi arcs. Then, one can provide the following ansatz for the "renor- 
malized" electron Green's function [14]: 

G{k,u)^ \ (C.23) 

where 

V{k) = U + 2t{cos + cos ky). (C.24) 

The first term comes from the fiuctuations of ag which is approximated by 
an on-site attraction: 5H — —Ud^d — —U/2{il)^h){K^ilj). The second term 
is due to the fuctuations of \Xij\ which induces 5H — —2tdjdi (this repro- 
duces the original hopping). The value of U then is fixed through the LSR 
/q°° o?a;/(27r) / d^k/ {2n)^'ismG — 5. Within this scenario, the antinodal direc- 
tions survives the binding energy V{k). On the other hand V{k) produces a 
negative chemical potential resulting in small hole pockets [14]. 

A recent progress in this direction has been done in Refs. [204,205]. 

It should be noted that the nodal-antinodal dichotomy emerging from the 
SU(2) save-boson analysis is distinct from that arising in the weak-C/ limit from 
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the RG approach; in the latter case, the antinodal regions are subject to promi- 
nent Umklapp scatterings that totally suppress the charge compressibility in 
the antinodal regions and produce a 2-gap scenario for underdoped cuprates 
(see Sec. 4). Additionally, in the latter scenario, the "pairing glue" around the 
nodal directions would stem from an Andreev-type scattering mechanism gen- 
erated by the hot (antinodal) regions [45,82,94]. Finally, one must bear in mind 
that within the weak-coupling RG analysis, the appearance of a truncated 
Fermi surface strongly depends on the form of the tight-binding Hamiltonian 
(one requires a next-nearest neighbor hopping) while details of the electron 
band structure are rather unimportant for large interactions. 

C.4 Theory of Quasiparticles in the superconducting state 

Finally, wc find it useful to discuss the superfluid density in the d-wavc su- 
perconducting state, following Lee and Wen [102]. We will also address the 
theoretical predictions from the SU{2) gauge theory [206]. 

Let us start with a superconductor where the elementary excitations in the 
superconducting state are well defined quasiparticles with dispersion: 

^(^) = [(^^-/^r + A|]^/^ (C.25) 
where represents the band dispersion and the d-wave gap is defined as: 

^ {cos - cos ky). (C.26) 

Here, following, Lee and Wen [102], we assume that there is a single gap 
scenario which is maximum of Aq at (0,7r). The important point is that in 
the d-wave superconducting state, there are four nodal points. In the vicinity 
of the node near (7r/2,7r/2), one has the anisotropic Dirac spectrum: 

E{k) = {vlkl + vlkiy/^, (C.27) 

where ki — {kx + ky — 7r)/V^, k2 — {kx — ky)/V2, v± oc Aq, and vp is the Fermi 
velocity. The presence of a vector potential shifts the quasiparticle spectrum: 

E{k, A) = E{k) - ^j{k) ■ A, (C.28) 
where j{k) is the current carried by the normal state quasiparticle with mo- 

— * — * 

mentum k. This equation is correct to first order in the vector potential A. In 
Rcf. [102], j is assumed to be —evp = — eS^e. The authors of Ref. [104] have 
introduced the Fermi liquid correction to the quasiparticle current [207]: 

j{k) = -eavp. (C.29) 
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a ~ 1 is related to the appropriate Landau parameter(s). The superfluid tensor 
is defined by: — — — ^A^. It can be decomposed as: 



pr = Ps(T = o)5^.-pr- (C.30) 

where ps{T = 0) is directly measured by the weight of the Drude peak in the 
normal phase and by the London penetration depth in the superconducting 
phase. Experiments [101] report that Ps{T — 0) oc S in agreement with the 
renormahzed mean-field theory as well as the slave-boson theory. On the other 

— * 

hand, p!^^ is given by the quasiparticle response to the A field. 
This can be evaluated by writing the free energy (at finite T), as 

F{A,T) = -keTY^H^ + e-^^(*^'^)), (C.31) 

k,s 

— * 

and differentiating twice with respect to A. Exploiting the definition of the 
current carried by the quasiparticles, we get: 

m ^dA^dA.dE' ^ ■ ' 

k 

f is the Fermi-Dirac distribution. Using Eq. (C.28) and (C.29), one gets: 

p,(T) _p,(T = 0) 



m m 

where [104]: 



aT, (C.33) 



2 In 2 2 'f^F 



a'^. (C.34) 

71 V± 

This is the result of Eq. (72). The interaction between quasiparticles has been 
neglected in (C.31); this is justified for small temperature and A because the 
density of states of quasiparticles vanishes linearly with energy, in contrast to 
the case of a Fermi liquid. For small doping S, the quasiparticle excitation is 
an effective way of destroying the superconducting state by driving ps to zero. 

An approximate expression of Tc then can be given by extrapolating ps — 0. 
Assuming that the gap or v± is independent of 5 for small S, then one predicts 
that Tc follows Ps{T = 0) cc S. This agres with the Uemura et al. plot [101]. 
It is relevant to observe that this scenario for Tc oc 5 is independent of the 
mechanism that produces superconductivity (Bose condensation, Kosterlitz- 
Thouless transition,...); it just relies on the existence of a gap with d-wave 
symmetry in the superconducting state, thus quasiparticles proliferate at finite 
temperature. The proximity to the Mott state appears through ps{T = 0) oc 5. 

In the slave-boson approach, the electron is decomposed into a fermion and a 
boson. The charge is carried by 5 bosons. The difficulty is that at the mean- 
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field level, the superconducting state is described by a condensation of bosons 
and the quasiparticle dispersion is given by the fermion (spinon) dispersion. 
Since the vector potential A couples directly only to the bosons, the shift in 
the quasiparticle spectrum is reduced and a ^ 1. In the original formulation, 
a — S. Within the SU(2) slave-boson approach, by invoking the spin-charge 
recombination mechanism (through the ladder diagram). Wen and Lee have 
shown that it is maybe possible to find ps(T = 0) ~ 5 whereas « ~ 1 is large. 
(This experimental fact may also be consistent with the fact that the Fermi 
arcs only become superconducting below Tg in the underdoped regime.) 



D Optical Lattices 

D.l Fermionic Atom 

It should be noted that the ground state electronic structure of alkali atoms is 
relatively simple: all electrons but one occupy closed shells, and the remaining 
one is in an s orbital in a higher shell. The electrons have no angular momen- 
tum, and the main coupling between the nuclear spin and the electron in the 
absence of magnetic field is the hyperfine interaction. The quantum number 
for the total spin F = 7 ± 1/2, where I represents the nuclear spin, behaves 
according to the usual rules for addition of angular momentum. The splitting 
between the levels F = J -|- 1/2 and F = / — 1/2 is given by 5f = {I + l/2)w 
where w is defined from the hyperfine interaction Hhf = wI.S. For example, 
has a nuclear spin / = 4; it is instructive to observe that the experiments 
performed in Refs. [57,58] have been done with two magnetic sub-levels of the 
F — 9/2 hyperfine manifold. In the case of ^Li, the nuclear spin obeys I — 1 
and therefore the total spin can take the two allowed values F = 1/2 and 
F = 3/2; experiments can be, e.g., realized in the two lowest hyperfine states. 



D.2 Standing wave 

In this Appendix, we briefiy remind how two counterpropagating laser beams 
produce a ID standing wave. First, when an atom is subject to an electric field 
E{x, t) this gives rise to a dipole type Hamiltonian Hdip — —jl.E{x, t) -|-h.c.; /i 
is the dipole operator of the atom. The electric dipole moment is proportional 

— * 

to the electric field: (jl) ~ aE. The energy of an atom in an electric field may 
be calculated using perturbation theory and a change in energy due to the 
electric field is ~ — alFp. The motion of an atom is thus governed by the 
Hamiltonian H — /2m + V{x) and the pseudo-potential V{x) oc \E{x,t)\^. 
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It is interesting to remind that the polarizabihties of alkah atoms are larger 
than that of hydrogen by factors with range from 30 to 90. 

A ID standing wave can be produced out of two running counter-propagating 

waves with the same intensity, as follows. Let (xi, 2:2, 2:3) be three unit vec- 
tors pointing along x, y, z direction, respectively. The electric field of the two 
waves is given by, say, Ei oc e'^'^^x^ and E2 oc e^'^'^^x^. The sum of the two 
electric fields is Ei + E2 oc cos{kx)x3, and this produces a standing wave in x 
direction; this creates an optical potential V{x) oc cos^lkx) in one dimension 
with periodicity A/2, where A = 27r/k is the wave length. Using one (two) 
further pair(s) of laser beams propagating in y {y and z) direction(s), a 2D 
(3D) periodic trapping then can be realized. The depth of the lattice in each 
direction is controlled by the laser intensity. 

D.3 XX Z model for the plaquette system at 3/8 filling 

Let us consider two plaquettes at sites R and R + 1. The effective XY term 
is given by exchanging |4/j,2/j+i) ^ |2^,4^_,,i) throough tunneling to the 
virtual state |3ij, 3i^+i); we have used the notations |2) = 1^2) and |4) = |^'4). 
This results in the XY term: the prefactor for this process is proportional to 
gH'^/Ab and the parameter g is essentially independent of interactions. In the 
dilute limit — with either a few 1^2) or 1^4) states — the dominant term 
is this XY process resulting in d-wave superfluidity. Close to the 3/8 filling, 
one also gets an Ising exchange term. However, the prefactor is not exactly 
equal to gH'"^/ Ab. For example, one also has to take into account ferromagnetic 
interactions of the form |4r4r+i) |4r4r+i) mediated by the |5r, 3r+i) state, 
and |2_R, 2_R+i) ^ |2r, 2r+i) mediated by |3r, Ir+i). Those processes result in 
a non-monotonous behavior of the coupling g^ [62], as a function U. For quite 
weak interactions, one expects gz < (?. At a general level, d-wave superfiuidity 
requires that A;, > and that g gz] this implies <^U < 2.7t [62]. 
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